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Abstract 


Riemann’s zeta function is one of the most important and fascinating functions of Euler’s math- 
ematics. By analyzing the material of Riemann’s hypothesis, we divide our analysis into the 
function ¢(z) and the proof of the conjecture, which has very important consequences in terms 
of the distribution of prime numbers. The proof of Riemann’s hypothesis follows from the simple 
logic,that when two properties are related, i.e. these equations are zero i.e. ¢(z) = C(1— z) =0 
while they have the proven 1—1 property of Riemann’s function ¢(z), then the Hypothesis will be 
true, because it has been proved from before simultaneously that Riemann’s functional equations 
hold. Therefore, there will be apply ¢(z) = 1/2 {because ¢(z) = ¢(1 — z) = 0 and also because 
¢(z) as well as (1 — z) are 1-1}. This, being true, will give the direction of all non-trivial roots 
to all be on the critical line, with a value of 1/2 on the real axis. Furthermore, to strengthen 
the proof we consider ”5 important requirements” that clearly prove that 100% of the roots of 
¢(z) = 0 are only on the critical line. There is of course the case ¢(q* z) =0 with q > 1/2 in R, 
which has different critical lines, but all those who deal with the Hypothesis of Riemann, do not 
mention it, because they claim that this form was not defined by Riemann from the beginning. 
But if we generalize the hypothesis to ¢(q* z) = 0,q > 1/2 in R then we have infinite critical 
lines that originate from and are related to the baseline critical line of Re(z) = 1/2, but are 
not identical to it and lie within the interval (0,1). In summary and by analytically examining 
and other similar functions to ¢(z), we get zeros with Re(z) in the interval (0,1) for the cases 
C(q * z) = 0 where q > 1/2 in R that mentioned, for Generalized case G — ¢(z,q) = 0 with 
at least 2 general cases out of 5 which is the generalized sum including ¢(z) and end for the 
Davenport Heilbronn case, which has zeros of a periodic Dirichlet series with different critical 
lines, within the interval (0,1). We are likely to have other constructed cases with relevance of 
¢(z). Therefore, whether or not the Riemann Hypothesis, holds depends on whether or not we 
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accept the fact that the Riemann Hypothesis is or is not fully mentioned, to functional equation 


C(z) = 0. 


Part I. The function ¢(z) is 1-1 on the critical strip 


1.Theorem 


The Hypothesis of Riemann focuses on the point where we must prove that if s = Re(s)+ 
+Im(s)*I, 


I) The functions ¢(s) and ¢(1 — s) are 1-1 on the critical strip. 
II) The common roots of the equations ¢(s)— ¢(1 — s) = 0 they have Re(s) = Re(1 — s) = 1/2 
within the interval (0,1) and determine unique position, which is called critical line. 


I.1. Proof: 
The functions ¢(s) and ¢(1 — s) are 1-1 on the critical strip. For this we need to analyze when and where 
the exponential function n? is 1-1 when n € Zt. 


The exponential function n*: 

I) is 1-1 in each of these strips defined by the intervals 27k/In(n) and 27(k + 1)/In(n) “where k € Z. 
For this we must prove two cases. If AC CA f: A — C then if f(z1) = f(ze) then z, = 22 that is.. 
t= 222A y1 = yo. Indeed, if 2] = 41 +2: Y1, 22 = 2 +7- yo are two points within into a such strip such 
that e*! = e*2 then, 


ntityrt — pxityit \n**| |n7? | are | = |n22n%*| >7,=2% 


and since #1 = #2, the relation e71™(™ = €7219(™) gives e* 4M) = er 92M), 50 y, — yo = 2kr/In(n), is an 
integer multiple of 27. But z; and z2 belong to the strip, so |y1 — y2| < 2a7/In(n). That is, the difference 
Yi — Yo is at the same time a multiple of 27/In(m) and at an absolute value of less than 27/In(n). The 
only case that this is true is when y; = yz. We finally conclude that z; = z2, so n? is 1-1 in the strip 
{z : 2ka/In(n) < Im(z) < 2(k + 1)7/In(n)}. ”lines down closed — open up”. We also notice that n* is on 
C — {0}. Because if w 4 0 and we put z = In|w|/In(n) + cargw/ n(n) then, 


nz= pln |wl+iargw/ In(n) _ nt lel pierge/ In(n) _ |w| ‘ niarge/ In(n) _ w 


2(k + 1)zi/tn(n) 


2kxi_ /tn(n) 


Figure 1: The z, and z2 belong to the lane width 27/ n(n) 


Formation of the strips 1-1 for ¢(s) and ¢(1 — s). 


al.) If we accept the non-trivial zeroes on critical strip of the Riemann Zeta Function ¢(s) as 81 = 01 + ity 
and s2 = 02 + it with |s2| > |si|, and if we suppose that the real coordinates 01, o2 of each non-trivial 
zero of the Riemann Zeta, [{1,2,3,4] function ¢(s) correspond to two imaginary coordinates t and ft, then, 
we have the following equations group: 


; 1 1 1 1 

C(o1 aes ti) _ poiteti 20 Joiti-ti : Bortiti END noitets eae 0 
. 1 1 1 1 

C(o2 aly te te) _ 1] c2tt-te Jo2ti-te : 302+1-te ae nr2tt te aa 0 


Taking the first equation and deducting the second, we obtain: 


[oe) 


: F 1 1 
G(r 7 ti) C(o2 Pas ta) = oy [os =a) . 
n=1 


co ne +t-t2 __ no +i-ty co n°? ev ta In(n) a= not ert In(n) ‘ 
~ nr2ttte . q_oiti-ti — Ss nratite ._oitity = 
n=1 n=1 


From the previous relation, we conclude that if 


; ; 2kr 
n7tita ptt 40 n>1 then og =o, Ate=t,+ 


int)’ S15 


That is to say t; and tz can take any value, but according to the previous relation. So from (Fig.1), that 
means ¢(s) is 1-1 on the lane of critical Strip in Intervals such as defined, therefore and on the critical Line. 


So n~* is 1-1 in the strip and {z : 2km/In(n) < Im(z) < 2(k + 1)7/In(n)}.” lines down closed — open up”. 


a2.) If we do the same work with |s2| > |s1| for the case ¢(1 — s) will we have: 


. 1 1 
c( oy ia ti) c(1 02 es ta) _ » (sea aan) ~ 


n=1 


co ni-72—tte ni-a—ttr ni-e2e tte In(n) a] ni-1 evr In(n) 


Co 
» 1—o2—it I—-o1-ity 1—oo—tt Toit mete) 
n 2 Fas n 1 A: n 2 Yao n 1 1 


n=1 n=1 


From the previous relation, we see that if 
nin g2—tte : Pe # 0,n > 1 
2k 
l-—o,=1- a A tg =t) + —,k= (1,2,... 
&({1— 02 a}Voz=a)An=h In(n)’ (L225) 


So, from (Fig.1), it follows that ¢(1 — s) is 1-1 on the critical Strip, So n~“~*) is 1-1 in the strip and 
{z: 2ka/In(n) < Im(z) < 2(k + 1)7/In(n)}.” lines down closed — open up” 


As we have seen before, the complex exponential form n* but and n~* is also periodic with a period 


of 27/In(n). For example n7+?°7* = n? -n*7*/lo8(") — ny is repeated in all the horizontal strips with 
2n/In(n) on (Fig.2) below. 
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Figure 2: All the ”bottom - closed upper “horizontal strips with 27/In(n) for n~*, n> 1, n € Z. 


nm 1 
If we are asked to find the strips that are 1-1 of the function ¢(z) = }) — it will be the union of the 
i=1 2 
strips formed with period 27k/In(n) where k,n € Z. This results from the analysis we made for the cases 
al, a2. Therefore, if we analytically assume that n = 2,3,... then the exponential complex function will be 


1-1 in each of the strips, ie. in each one set of the form le 


Sn = {z: 2kr/In(n) < Im(z) < 2(k + 1)n/In(n)},k € Z 


co 
and therefore in each one subset of the unified expanded set S = J S;. The largest strip obviously contains 
i=2 


all the rest it is wide, and has width 27/1In2 i.e. with n = 2. 


I.2. Overlapping contiguous Strips. 

If we start from the upper bound of the strip with n = 2 and k = 1, Fig.3 the upper ones will descend to 
zero by increasing n with a width difference the interval 6 = 2a/In(n — 1) — 27/In(n). In addition, we will 
prove that the strip i = n > 2 primary (k=1) that generated in ascending order, the upper limit of that 
created of its double width enters the zone of the previous of (i = n —1,k = 1) also primary. With this 
logic we can find for which n we have overlapping strips in relation to the increase of n on the imaginary 
axis. Below we find for which n this dimension applies which is useful to distinguish them but also to make 
them visible schematically. 
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Figure 3: The upper ones will go down to zero by increasing n 
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talked about integers then n > 3. The same phenomenon is also created at higher level strips. 


and because we 


Apply 27/In(n —1) < 2-:2-7/Inn>n< (n-1)? > n?-3-n+1>053n> 


“1 
1.3. The) ¢(z) = S- — | it is directly 1-1 from the comparison of 2 points 2, z2. 
i=1 


Proof: 
From the equality 


1 tl 
— = SS HC (zq) 1 FQ) 4 grt) 4 pn tm) = ¢(z9) = 


22 


3 


= 1 4 2Q- (ety) 4 g-(@wetyet) 4 4 (watyet) 


with comparison of similar terms the result arises. 


2:-7-k 
1 =a. Ay, = y2 - ——— kez 


In(n) ’ 


After to the analysis of cases before, will be arises 1-1 as shown in sections {I,.a1, I1.a2} pages 1& 2. Also 
according to Lagrange’s generalized theory, each of the roots according to the obvious relation ¢(z) = ¢(1—z) 
will result from the generalized theorem of Lagrange. From the equation that results from 3 Functional 


Equations, with condition of common roots. We take ¢(z) = ¢(1 — z) and for first approach we then follow 
the 3 first terms on each ¢-equation we will have in the Regular form ¢(z), ¢(1 — z) with the analysis below. 


1 1 
1-2 42-%43-%=1 qa z49 qa 243 qa 239 743 Sao tae 


3 oY 
yaaa ee ae ee ae Ses ee =O Oth ate 3 (3) =0 
With the replacement 


Log(y) +2-k-+a-i 
Log(6) 


67 =y => z- Log(6) = Log(y)+2-k-t-isz= 


Log(y)+2-k+a-4 
Zin = 
Log(6) 


we consider it as initial value for the solution of Transcendental equation with the method of Lagrange 
inversion, then: 


w 


Log(y) +2: kn sans ALLA, 


_ — Logly) +2 km a , SS (-1)"— avo! 1 [3.4 Log) _9.f(2 Log(6) 
tos) * Tw ye (sea) ie : (3) we 


ys 


Some solutions for k € Z (Table 1) 


18t 0.4985184869351 -/+  3.6072014300530 I, 
2°¢ (.4997077392995 -/+  10.329023684711 I, 
3'4 0.5003426560276 -/+  13.990255078243 I, 
-/4 
-/4 


4‘ (.4993087825728 + =: 17.585326315022 I, 
5' 0.5026747758005 + =. 21.304192633721 I, 


TTD 
Oo» Wr 


I 


Table 1: 6 Roots of equation ¢(z) = ¢(1 — z) 


And so on for the sequence of infinity. Now we see that we approach the real roots and its complexities. 
With the additional method Newton or the Bisection method we approach more rapidly to them. And these 
roots will be contained in corresponding strips that we have previously defined and in closer relation with 
their position on the critical line. 

By this logic for the previous approximation equation in relation to its solutions of, ¢(z) = 0 the strips 
which determine the imaginary roots of the equation, will be between the intervals 2-7-k/In6, k € Z as 
defined by the analysis of the generalized Theorem of Lagrange. The strips of the imaginary part of the 
roots, as they appear for 6 consecutive intervals of the approximate equation with method of Lagrange that 
we have solved above and we will have the results on the (Table 2). 


Strips — Imaginary Roots: 
18* (7.013424977055179 — 3.506712488527589) 
24 (10.52013746558276 7.013424977055179 
3™¢ (14.02684995411035 10.52013746558276 
4*® (17.53356244263794  14.02684995411035 
5 (21.04027493116553 —17.53356244263794 
6 (24.54698741969312  21.04027493116553 


OS EE) 


Table 2: The strips of the imaginary part of the roots 


It is clear the first root 3, 607..i is located in 1** of strip (7.01 — 3.506), the second 10.3290..i in the 2** strip 
(7.013 — 10.52) etc. Now the root 13.99..2 which approximates the first root of ¢(z) = ¢(1— z) =Oie. the 


14.1347..i, is in the 4%* strip (14.02 — 17.533) with a lower limit the 14.026..i a value that is very close to 
the one required for the approximate and exclusive (1** root) when ¢(z) = 0 of the Riemann Hypothesis. 


Corollary 1. 
The only roots of the Zeta function not included in the set {z € C : 0 < Re(z) < 1} are the points 
—2,—4,—6.. [11. p47] 


Proof 

On the functional equation: 

¢(1 — z) = 2-¢(z)- T(z) - cos(a/2- z) - (27)~* we know that for Re(z) > 1 the functions ¢(z) and T'(z) do 
not equal zero (Proof-Th.3.p6). We also find that for Re(z) > 1 the 1 — Re(z) < 0 and putting u = 1- z 
we will find all the roots of ¢(u) for Re(wu) < 0. Therefore we will have that ¢(w) will be zero, where it is 
zeroed the i.e. z = 3,5,7. Then the roots of ¢(w) for Re(u) < 0 will be the points —2, —4, —6... and will be 
all the roots of the function ¢ out of the strip {z € C: 0 < Re(z) < 1}. 


Theorem-Helpful 1 
The Riemann’s Z-function has no roots on the lines Re(z) = 1 and Re(z) = 0. 
The proof is detailed in the book [11. p50-51). 


II. The common roots of the equations ¢(s) — ¢(1 — s) = 0 they have Re(s) = Re(1 — s) =1/2 
within the interval (0,1) if moreover apply ¢(s) = 0. (Refer to Theorem 3, p.8) 


Proof: Let us assume z to be such that for complex z = xp + iyo, 0 < Re(z) < 1A Im(z) 4 0 and 
¢(z) = ¢(1— z) = 0. According to the two equations, they must apply to both, because they are equal to 
zero that: ¢(z) = ¢(Z) = OAC(1— z) = ¢(1— z) = 0. But from [Theorem 1, II, al, a2] the ¢(z) and ¢(1— z) 
are 1-1 on the critical Strip. If suppose generally that ¢(ap + yo -2) = ¢(1— 29 £ yf - 2) = 0 then: 


Yo = Yop AZo = 1-2, 


but because are the 1-1 then we will apply two cases for complex z i.e. 


= 


Nl rR 


/ / / / Se ae 
Ly =1—Xo, LAX, YH Yo VX =1—2o0, LT =X = 


{Corollary 1, Theorem Helpful 1}, because we suppose ¢(z) = ¢1— z = 0.] We conclude that 0 < x < 1. 
Therefore we have three cases: 


1 
IL.1.. If | 2.25 = 5 


Ayo = yo. Which fully meets the requirements of hypothesis! 


we apply the obvious i.e. ¢ (1/24 yo-t) =¢(1-—1/2£y-t) =¢ (1/24 y- i) =OA 


1 
II.2. If 0<2#mA#A5<1 in this case if zo < x Aa +29 = 1 or 2p > 2 AX +25 = 1 then for the 


Functional equation ¢(s) — ¢(1 — s) = 0 simultaneously apply: 
C (xo + yo 4) = ¢ (to — yo: t) =ONC (x + YO ° 4) =¢ (29 — Yot) = 0. 
So let’s assume that: 
1 no <a Am +a =14 m= 5-ahm=5+daxb, But apply v9 +29 =1>1+b-a=1=> 


1 
b=a. That is, a,b is symmetrical about of 5° But the ¢(z) and ¢(1 — z) are 1 - 1 on the critical Strip 
then apply the relation: 


1 ak 1 1 
=-—aty-i} =C(=+aty-i}) Ayo=yo > must --a==+a>2a=05a=0 
2 2 a : 2 2 


therefore xo = x 
1 1 

24 go > th At0+2=1> 2= 5 +aA% = 5 a as before. Then because the ¢(z) and ¢(1 — z) are 

1-1 on the critical Strip will apply: 


1 ; 1 Pas ! 
G 3+ a=Yo% =C(5—4+%-2 A Yo = Yo 


1 1 
that -+a=--—a=>2a=0=>a=0 therefore r 9 = xp. We see that in all acceptable cases is true that 


” The non-trivial zeros of ¢(z) — ¢(1 — z) = 0 have real part which is equal to 1/2 within the interval (0,1). 
In the end, we see still one partial case that appears in the roots of the equation ¢(z) —¢(1 — z) = 0. 


11.3. If] 1 <2 Ax < 0| furthermore, in this case with that is, it is symmetrical about rp + 7) = 1 then 


for the functional equations we apply: 
¢ (0 — yo: 4) = ¢ (1-20 — yo 4) AC (xo + yo: t) = ¢ (1 — 20 + Yo) A Yo = Yo 


and for the three cases it is valid yo = 0 V yo # 0. This case applies only when ¢(z) = ¢(1 — z) OA ao = 
= 1+4a, x) = —a, and not for the equality of ¢ function with zero. But in order to have common roots of 
the two functions ¢(z) and ¢(1 — z), «because they are equal to zero and equal to each others, then first of 
all, it must be i.e 29 = xj = 1 — 2 and also because the ¢(z) is 1-1 and for each case applies (yo = 0 and 


yo = 0) or (yo #0 and yj # 0), therefore we will apply xo = x = 5 A Yo = Yo and with regard to the three 


cases {II.1, II.2, I1.3} the cases {II.2, II.3} cannot happen, because they will have to be within (0,1) and 
will therefore be rejected. 


Therefore if ¢(z) = ¢(1— z) = 0 then because ”these two equations of zeta function are 1-1 on the lane of 
critical strip” as shown in theorem 1, is “sufficient condition” that all non - trivial zeros are on the critical 
line Re(z) = 1/2. That means that the real part of z of ¢(z) = 0, equals to 1/2”. So Theorem 1 & Corollary 
1, they have been proved and also helps Theorem 3 below. 


2. Theorem 2 
For the non- trivial zeroes of the Riemann Zeta Function ¢(s) apply 


There exists an upper-lower bound of Re(s) of the Riemann Zeta Function ¢(s) and more specifically 
In2 Inz 
In2n’ In2r 
lower bound are distributed symmetrically on the straight line Re(s) = 1/2. The average value of the upper 
lower bound of Re(s) = 1/2. 


in the closed interval | | . The non-trivial zeroes of the Riemann Zeta function ¢(s) of the upper- 


Proof: 
Here, we formulate two of the functional equations from E - q. Set 


¢(1 — s)/¢(s) = 2(27)~* cos(7s/2)T(s), Re(s) >0 
¢(s)/C(1 — s) = 2(27)*—! sin(ws/2)T(1—s), Re(s) <1 
We look at each one equation individually in order to identify the set of values that we want each time. 


2.1. For the first equation and for real values with Re(s) > 0 and by taking the logarithm of two sides of 
the equation [5], we have: 


C(L— s)/¢(s) = 2(2m)~* cos(/2)P(s) + 


log[¢(1 — s)/¢(s)] = log[2] — s log[27] + log[cos(as/2)T'(s)] + 2kxi 


but solving for s and if f(s) = cos(ms/2)T(s) and from Lemma 2 [6, p.6]: 
if lim ¢(1 — s)/¢(s) = 1 or log lim ¢(1 — s)/¢(s)| = 0 we get s = 
S—> So S780 


log|f(s)] + 2kri 
He ( log[27] 


log[2] _ log[f(s)] + 2kmi 
log[27] ° log[27] 
log[2] 
log/27] 


with 


= 0.3771 |. This is 


) > 0, finally because we need real s we will have | Re(s) > 


the lower bound, which gives us the first ¢(s) of Riemann’s Zeta Function. 


2.2. For the second equation, for real values if Re(s) < 1 and by taking the logarithm of the two parts of 
the equation, we will have: 
¢(s)/C(1 — s) = 2(27)*“" sin(as/2)T(1 — s) > 
log[¢(s)/¢(1 — s)] = log[2] + (s — 1) log[2z] 4+ log|sin(a - s/2)[(1 — s)] + 2kai 


but solving for s and if f(s) = sin(ws/2)T(1 — s) and from Lemma 2[6, p.6] if lim [¢(s)/¢(1 — s)] = 1 or 


log|[t] — log[f(s)| +2kai log[f(s)] + 2kmi\ — _ 
log[27] log[27] Withee ( log[27] ) are 


log jim ¢(s)/¢(1 — s)| =0 we get s = 


1 
In the following, because we need real s we will take. | Re(s) < cael 
og [27 


= 0.6228 |. This is the upper bound, 


which gives us the second of Riemann’s Zeta Function of ¢(s). So, we see that the lower and the upper 
bound exist for Re(s) and they are well defined. 


logiz] 


. : : . I 
2.3. Assuming that sz = Glow + it, and s, = upper + it, with o1low = log]27] 


Oo 
log[27] 
apply 01 = Olow and 02 = Gupper , (Fig.4). If we evaluate the difference and 


and Oupper = we 


Re Al, = Re s', — 1/2 = oupper — 1/2 = 0.1228 
Re As, =1/2— Re s, = 1/2 — oipy = 0.1228. 


This suggests for our absolute symmetry of 


it Oy tity 


Figure 4: Arrangement of low and upper of real part of zeros from 2 functional equations of Riemann. 


2.4. The average value of the upper lower bound is Re(s) = 1/2 because from (Figure 4): that is ap- 


ply 
log[2] + log[z] 


Re(s) =1/2- log[2m] 


=1/2 


3. Theorem 3 
The Riemann Hypothesis states that all the non-trivial zeros of ¢(z) have real part equal to 1/2. 


Proof: In any case, we Assume that: The Constant Hypothesis ¢(z) = 0, z € C. In this case, we use 
the two equations of the Riemann zeta function, so if they apply what they represent the ¢(z) and ¢(1 — z) 
to equality. Before developing the method, we make the three following assumptions: 


3.1. Analysis of specific parts of transcendental equations, which are detailed 


3.1.1. For z € C, {a* =0=>a=0/A Rez > 0} which refers to the inherent function similar to two Rie- 
mann zeta functions as (27)~* = 0 or (27)~+ = 0 and it seems that they do not have roots in C—Z, because 


(27) £0. 
3.1.2. This forms ['(z) = 0 or ['(1 — z) = 0 do not have roots in C — Z. 
3.1.3. Solution of sin(z/2z) = 0 or cos(a/2z) = 0. More specifically if z = x + yt, then. 


1. sin(a/2z) = 0 = cosh[(z - y)/2] sin[(a - #)/2] + I - cos[(a - x) /2] sinh[(z - y)/2] = 0 
If k,m € N (Integers) then all the solutions can be found with a program by the language of Mathematica. 


{a = (2+4-k),y=2-i-(1+2-m} 
{r=4-k,y=2-i- oe m} 


i) 
ii) 
ii) {a =(-14+4-k),y=i-(-1+4-m} 
w) {a =(-14+4-k),y=7-(1+4-m} 
v) {©=(14+4-k),y=1- ee 
vi) {#=(14+4-k),y=i-(1+4-m} 
vit) {w= (24+4-k), y=i- (4-m} 
viti) {a =(24+4-k),y=2-i-(14+2-m} 


From the generalized solution, it seems that in the pairs (2, y) will always arise integers which make impos- 
sible the case 0 < x < 1, therefore there are no roots of the equation sin(7z/2) = 0 in C — Z. 

Also, because as we see, all the roots given by the union of the sets x > 1V x2 =-—1Va=0 all are Integers. 
So it cannot be true that for our case: 


cos(7/2z) = 0 = cosh|(ry)/2] cos[(aax)/2] — I sin[(m - x) /2] sinh[(ay) /2] = 0 


If k,m € N (Integers) and the solutions are: 


{x =4-k,y =i-(-1+4-m)} 


) 

) 

:) 

iw) {es (-144-Biy = 

v) GS ay Saeko 

vi) {#=(1+4-k),y=i-(1+4-m} 
vit) {v= (2+4-k),y=i-(-1+4-m} 
viti) {a =(2+4-k),y=i-(14+4-m} 


As we can see, again from the generalized solution, it seems that in the pairs (x, y) will always arise inte- 
gers, which make impossible the case is 0 < x < 1, therefore there are not roots receivers of the equation 
cos(7z/2) = 0. Since all the roots are given by the union of the sets x > 1 V2 = —1V a =0 it follows that 
they are Integers. 


Great Result 

<With this three - cases analysis, we have proved that the real part of x, of the complex z = x + yi cannot 
be in the interval 0 < x < 1, and in particular in cases {3.1.1, 3.1.2, 3.1.3}, when they are zeroed. Therefore 
they cannot represent roots in the critical line>. 


3.2. Therefore, now we will analyze the two equations of the Riemann zeta function and we will try 
to find any common solutions. 


3.2.1. For the first equation and for real values with Re(z) > 0 we apply: 
Where f(z) = 2(27)~? cos(m - z/2)I'(z) > C(1—z) = f(z)-¢(z) but this means that the following two cases 
occur: 


3.2.1.1. If] ¢(1 — z) = ¢(z) |, where z is complex number. 

This assumption implies that ¢(z) - (1 — f(z)) = 0 => ¢(z) = 0. In theorem 1.(I , II), (pages 1 — 8), 
we showed that the functions ¢(z) and ¢(1 — z) are 1-1 and therefore if ¢ (ap + yot) = ¢ (ap + yp2) then, 
but we also apply that ¢ (ao + yot) = ¢ (ao — yot) = 0, because we apply it for complex roots. The form 
¢(1 — z) = ¢(z) means that if z = xo + yo i we apply 1 — x = vo namely xp = 1/2 because in this case 
it will be verified that ¢ (1 — xo — yt) = ¢ (1/2 — yt) = ¢ (1/2 + yot) = ¢ (1 — 20 + yOt) = 0 and because 
Yo = Yo = C (x0 — yot) = ¢ (1/2 — yot) = ¢ (1/2 + yot) = ¢ (ao + yot) = 0 and these are the two forms of 
the ¢ equation, they have common roots and therefore it can be verified by the definition of any complex 
equation, when it is equal to zero. Therefore if z = x) + yo i then xp = 1/2 which verifies the equation 


¢(z) =0. 


3.2.1.2. If} ¢(1 — z) # ¢(z) |, when z is a complex number. 
To verify this case must be: 


G(z) = C1 — 2): f(z) > C(z) = 0A f(z) =0. 


But this case is not possible, because as we have shown in Section 3.1 (3.1.1, 3.1.2, 3.1.3), the individual 
functions of f(z) cannot be zero when z is a complex number. 


3.2.2. For the second equation and for any real values with Re(z) < 1 we apply: 
¢(z)/C(1—z) = 2(27)*~? sin(-z/2)P(1—z) => C(z) = f(z)-¢(1—2), where f(z) = 2(27)*~1 sin(-z/2)T'(1—z) 
but this also means that two cases occur: 


3.2.2.1. If] ¢(1 — z) = ¢(z) |, when z is a complex number. 


This case is equivalent to 3.2.1.1, and therefore if z = ap + yo i then xp = 1/2 in order to verify the equation 
¢(z) = 0 we follow exactly the same process algebraically. 


3.2.2.2. If |¢(1— z) 4 ¢(z) |, where z is a complex number. Similarly, the above case is equivalent 


to 3.2.1.2 and therefore it cannot be happening, as it has been proved. 


4. Forms of the Riemann ¢ Functional Equations 
The Riemann ¢ function has three types of zeros: [7] 

Eq. Set ¢1—s) = 2-(2-2)~*-cos(r-s/2)-T(s) + ¢(s), Res > 0 

ms C(s) .(2-m)8-} . sin(a-8/2)-T(1—s)-¢(1—s), Res<1 


oe) 


(I, II) 


I 
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And 


1-—s 
ile 
_ s—1/2 2 
¢(s) = ne-¥/2. A 2 Les) (111) 
(3) 
2 
Usually referred to as the trivial zeros, and non-trivial complex zeros. Therefore has been proved that any 
non-trivial zero lies in the open strip {s € C : 0 < Re(s) < 1} that is called the critical Riemann 


Strip: And all complex zeros of the function ¢ lie in the line {s € C : Re(s) = 1/2} which is called the 
critical line. 


Solving of the Riemann ¢ Functional-Equations. 


To find the imaginary part we must solve the functional equations (Eq. set), and for the cases which 
the real part of the roots lies on the critical line Cases 4.1, 4.2 [6]. 


4.1. 18* type roots of the Riemann zeta functions (15 equation from the Eq. set) 
For the first category roots and by taking the logarithm of two sides of the equations, and thus we get [5]: 


(1 — z)/¢(z) = 2(27)~* cos(m - z/2)T(z) > 
=> log[¢(1 — z)/¢(z)] = log[2] — z log[27] + log[cos(a - z/2)T'(z)] + 2kai 


and the total form from the theory of Lagrange inversion theorem, [5] for the root is pi(z) = s which means 
that f(s) =p; '(s) = z, but with an initial value for s which is 


log[2] 2-k-m-t 
log[27r] log(27) 


Sin 


by setting values for k we can simply calculate the roots of the (18* equation from the Eq. set, (Table.3)). 


18t 0.377145562795 ~—-/4 

2°¢ 0.377145562795  -/+  3.41871903296 I 
3'¢ 0.3771455627955 -/+  6.83743806592 I 
4*® 0.3771455627955 -/+  10.2561570988 I 
5 0.3771455627955 -/+ 13.6748761318 I 
6 0.3771455627955 -/+ 17.0935951646 I 


Table 3: Solution of 1* equation from the Eq. set, The firsts 6 Roots. 


A simple program in mathematica is: 


k > 4; t = (log[2] + 2k7i) / log[2z]; 

z=s+ 3) ((—1/ log[2z])” /Gammal[w + 1]) + D [(s’) - log[cos|z - ¢/2]]+ 
+log[Gamma|s]] — log|Zetta[1 — s]/Zetal[s]])”, {s, w — 1}]; 

N{z/.s > t, 10] 


With k = 0, 


eR 
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But because the infinite sum approaching zero, theoretically s gets the initial value 


log[2] 


2-k-m-1t 


log[27r] 


log(27) 


So we have in this case, in part, the consecutive intervals with k = n and k = n-+1 for any n > 4 and for 


the imaginary roots. 
4.2. 24 type of roots of the Riemann zeta func 


Same as in the first category roots by taking the | 
get: 


tional equations (2"¢ equation from the Eq. set) 


ogarithm of two sides of the equations, and thus we 


¢(z)/C(1 — z) = 2(2n)** sin(ax - z/2)T (1 — z) > 
=> log[¢(z)/¢(1 — z)] = log[2] + (z — 1) log[2z] + log[sin(a - z/2)[(1 — z)] +2-k- 7-4 


Now, we will have for total roots of 3 groups fields (but I have interest for the first group), and therefore 


for our case we will get pi(z) = s which means that: 


f(s) =p 


but with an initial value for s that is 


1 (8) =z, 


log[7] 2-k-em-t 


in 


log[27] log (27) 


the overall form from the Lagrange inverse theory succeeds after replacing the above initial value and 
therefore for the first six roots after calculating them we quote the following (Table.4). 


15* 0622854437204 
2nd () 622854437204 
3° 0.622854437204 
4th (.622854437204 
5*h 0.622854437204 
6° 0.622854437204 


-/+ 
-/+ 3.41871903296 I 
-/+ 6.83743806592 I 
-/+ 10.2561570988 I 
-/+ 13.6748761318 I 
-/+ 17.0935951648 I 


Table 4: The firsts 6 Roots of 2"? equation from the Eq. set 


A simple program in mathematica is: 


k > 4; t = (log[x] + 2kni)/log[2m}; 


+log|Gamma|{1 — s]] — log[Zetta|s] 
N{z/.s > t, 10] 


z=s+ 3, ((—1/ log[27])” /Gammal[w + 1]) + D [(s’) - log[sin[z - s/2]]+ 


/Zetal1 — s]])”,{s,w — 1}]; 


Therefore, with k = 0,+1,+2,... 


But because the infinite sum approaches zero, theoretically s gets initial value 


log[7] 


2-k-m-t 


log[27r] 


log(27) 
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So we have in this case, in part, the consecutive intervals with k = n and k = n-+1 for any n > 4 and for 
the imaginary roots. And for the cases we have for Im/(s) the relationship 


_ 2-a-k 


We) ouome 


For the 3'¢ functional equation, one has been previously put on the (as in the other two in 4, page 10) and 
we take as imaginary part 


Following to complete the roots of the two sets of the functional equations Eq. Set, we solve the functions 
as cosine or sin according to its Generalized theorem of Lagrange [5]. 


4.3. Transcendental equations for zeros of the function (Explicit formula) 


The main new results presented in the next few sections are transcendental equations satisfied by indi- 
vidual zeros of some L-functions. For simplicity, we first consider the Riemann-function, which is the 
simplest Dirichlet L-function. 


sym otic equation satisne e n- zero on e critical line. 
* Asymptotic equati tisfied by the n-th the critical line. [9, 10, 11] 


As above, let us define the function x(z) = 1~*/?T'(z/2)¢(z) which satisfies the functional equation a(z) = 
a(1—z). Now consider the Stirling’s approximation I'(z) = /2mz77-1/?e-? (1+ 0 (z7')) where z= x +iy, 
which is valid for large y. Under this condition, we also have 


2 = exp (i(ylogy +) + «logy + +2+O(y")). 


Therefore, using the polar representation ¢ = |¢\e’®"8S and the above expansions, we can write as 


A(e,y) = vBmm2!2 (ZY 8" e-ev/Ate(e + iy (140 (2), 


a (e 1) + arg ¢(a + iy) +O (y+). 


ren) Bi (3h) + 


The final transactions we end up with 
y y eee csate. 1 
= + og (54) lim — ee 
27 oe 27e = + lim, © argc (5 + aa 
Establishing the convention that zeros are labeled by positive integers, z, = 1/2+7-y, where n = 1, 2,3,4,..., 
we must replace n + n—2. Therefore, the imaginary parts of these zeros satisfy the transcendental equation 


n 


t t 1 1 11 
Eq. i | +t ii ~4+6+ity) =n-—. 
oe 27 os (5) jim, Saree (5 + ee ) . 8 


Let us recall the definition used in, namely: 


1 1 1 1 
Stes tee 1 eer 1 : 
so ip benc(}-6) =a, 2 [ne ($604) 
These points are easy to find, since they do not depend on the uctuating S(y). 
We have: ; : 

¢ (5 7 iv) =¢ (5 +iy) Gy), Gly) = ©) 
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The Riemann-Siegel # function is defined by 


1 4 = 
d= arr ( +5) y log Vr. 
4 2 
Since the real and imaginary parts are not both zero, at ys? then G = 1, whereas at yy? then G = -1. 
Thus 


S E (5 + vs”) = 0 for J (vi?) =(n-—1)z, 


R E € Fe iw) =0 for 0 () =(n—1)r. 


they can be written in the form of Eq.b. 


27(n — 7/8) (-) 27(n — 3/8) 


Eq. b jy = Wlet(n—7/8))' ™ ~ Wien —3/8)] 


where above n = 1,2,... and the Labert - Function W denotes the principal branch W(0). The ys?) 
are actually the Gram points. From the previous relation, we can see that these points (Fig. 5) are ordered 
in a regular manner [8, 12] 


YP < yi? 


(+) 


< Yo ae 


< Y2 < yf < yf eet 


Figure 5: Gram points 


With this method, we are able to create intervals so as to approximate the correct values of imaginary 
part of non-trivial zeros. We did something similar in the cases 4.2. What it has left, is to find a method 
that approximates the values of imaginary part accurately. 
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Figure 6: The black dots represent the zeros of ¢(s) function including 
possible zeros which do not lie on the critical line. 


The zeros of the Riemann Zeta function. 


“PROGRAMMING” 
5.1. The M Function — Bisection Method > ” 
Knowing the time of the successive steps (k,k +1) of the relationship of imaginary parts Im(z) = aD y? 
1 


2-a-k 
FN) = alee 


¢(1/2+7-y) =0 using the Bisection Method. 


(Eq I, II, Ill, page 10) with k € N, we can calculate the roots by solving the equation 


Bisection is the division of a given curve, figure, or interval into two equal parts (halves). A simple 


bisection procedure for iteratively converging on a solution, which is known to lie inside some interval [a, b] 


a+b Socal 


proceeds by evaluating the function in question at the midpoint of the original interval y = 


testing to see in which of the subintervals [a, (a + b)/2] or [(a + 6)/2], 6] the solution lies. The procedure is 
then repeated with the new interval as often as needed to locate the solution to the desired accuracy. Let 
dn, bn be the endpoints at the n'” iteration (with a, = a and b,; = b ) and let r, be the n™ approximate 
solution. Then, the number of iterations required to obtain an error smaller than ¢ is found by noting 


— it 
that... bn — Gn = ie and that r, is defined by rn = 5 (dn + b,). In order for the error to be smaller 


b 
Qn-1 


1 
than ¢, then... |rn —r| < 5 (by — Gn) = 2-"(b— a) < €. Taking the natural logarithm of both sides gives 
—n|n2 < Ine —In(b— a). Therefore, we have for steps: 


5.2. M-function of the Bisection Method.. 


We define the functions Mj, Mj on an interval (a,b) according to the scheme: 
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I. Mj, with Mq_, the Nearest larger of Mg_m where 


Met Men > 2,d>m>d-2, 
b 
M, =" | 


II. Mj, with Mg_1 the Nearest smaller of Mg_m where 


Mai * Mem 4 >2,d>m>d-2, 
b 
M, =" Ra 


For calculating the roots on solving the equation ¢(1/2+7-y) = 0 take the limit rg according to the scheme: 


rs = lim Mt = lim M,,s € N with Mf and M7, 


Qn-k Qn (k+1 
belong in the interval (a,b),d € N and also («= Tees alu N,k > 3 for the x of 


C(1/2+i-rs) =0. 


5.3 Program in Mathematica for the Bisection method of ¢(1/2+7-y) =0. 


2-7-k 2-m-(k+1 
Using the Intervals | a = a a pee tbetsl) and successive steps, we can compute all the roots 
log/27] log/27] 


of ¢(1/2 +7-y) =0. We can of course use three types such intervals more specifically in generall than 


2imek 2m (h ET) 


ag log[27]’ slog [27] 
_2emek 2-H (k+1) 
logit] 7 Mogi] 


3.a = yn b= yi. 


yO 27(n — 7/8) yO = 27(n — 3/8) 
- W [e-!(n — 7/8)]’ 7” W [e-!(n — 3/8)] 


We always prefer an interval that is shorter, in order to locate fewer non trivial zeros. The most important 
is to calculate all the roots in each successive interval and therefore only then we will have the program 
8-7 10-7 


= jog[2n]’ = jogi2m] and Error approximate tol = 107° and 


for data example: Integer k = 4, and (« 


Trials n = 22.. 


“Programm for Bisection method” 


A program relevant by dividing intervals.. 
This program gives very good values as an approach to the roots we ask if we know the interval. Selecting 


(« 2. ak a) 


= logi27]’ Ss dog [27] 


the interval for the case ,k € N,k > 3 the results are given below in 


(Table.5).. 
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Clear["**"]; 
[x_]:=Zeta[1/2+x*I]; 
=Input["Epilogh k"); 
a=(2*k)*z/Log[2a]; 
=2*(k+1)*2/Log[2]; 
o-Input["Enter tolerance"]; 
=Input["Enter total iteration"; 
=N[(-Log|tol]/Log[10]+4)|; 
f[Arg[f[c]]>0,{Print["No solution exists"] }]; 
rint("n a DC suse £(C)...00.""]5 
o[{c=N| (atb)/2,g]; 
f[Arg[f[c]]<0,a= c,b=c], Print[PaddedForm[i,10],PaddedForm[N[a],{7,7}], 
addedForm[N[b],{7,7}],PaddedForm[N|[c],{7,7}], 
addedForm|N|[f[c]].{7,7}]]If[Abs[a-b]<tol && Abs[N[f[c]]|<tol*1000, 
{Print["The solution is: ",N{c,g]] Exit{]}]},{i,1,n}]5 
rint["The maximum iteration failed,No solution exists"); 


1 
Calculation 18* Zetazero z, = 5 + 14.1347216.. 


b 
17.0935951600 
15.3842356500 
14.5295558900 
14.5295558900 
14.3158859500 
14.2090509800 


a 

13.6748761300 
13.6748761300 

13.6748761300 
14.1022160100 
14.1022160100 
14.1022160100 

14.1022160100 
14.1289247500 
14.1289247500 
14.1289247500 
14.1322633500 
14.1339326400 


c 
| 
13.6748761300 | 
15.3842356500 _| 
14.5295558900 | 
| 
| 
| 
| 


14.1022160100 
14.3158859500 
14.2090509800 
14.1556335000 14.1556335000 
14.1556335000 14.1289247500 
14.1339326400 
14 14.1343499700 14.1347672900 14.1343499700 
15 14.1345586300 14.1347672900 14.1345586300 | 
16 14.1346629600 14.1347672900 14.1346629600 | 
14.1347151300 14.1347672900 14.1347151300 
18 14.1347151300 14.1347412100 14.1347412100 
19} 14.1347151300 14.1347281700 | 14.1347281700 | 
20) 14.1347216500 14.1347281700 14.1347216500 | 


COIN IDM BLWIN [Re 


Table 5: The first non-trivial Root, Imaginary Part. 


With final value 14.1347216500, error near of 10~’. This value is the approximate root of the by the nearest 
error < 1077. 

In the event that we have two or more uncommon roots within the interval, we divide similar successive 
intervals in the order of finding of the first root either above or below. In such a case we have the k = 13, 
and at the interval (47.8620664, 51.280785) the two roots are 48.0051088 and 49.7738324. 


5.4 Explicit formula and the Zeros of ¢(1/2+7-y) =0. 
Consider its leading order approximation, or equivalently its average since < arg ¢(1/2+i-y) >> 0 [11, 12]. 


ta 
27e 


tn in 
Then we have the transcendental equation on log ( ) =n- re. Through the transformation t,, = 
7 
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11 
3 =): Comparing the previous results, 


we obtain 


11 
Qn (w = =) x, this equation can be written a 2,e™ = e+ (w “ 


where n = 1,2,3... 


5.5 Programm by Newton’s method, which finds the Zeros of ¢(1/2+7-y) =0 


Using Newton’s method we can reach the roots of the equation ¢(1/2 +%-y) = 0 at a very good ini- 
tial value from y > t, by the explicit formula. 


It follows a mathematica program for the first 50 roots by the Newton’s method (Table.6). This method 
determines and detects the roots at the same time in order to verify the relation ¢(z) = 0 and always 
according to the relation... 


where n = 1,2,3... and where W is the W-function. 


Following is the program in mathematica and the Table.6 of complex roots... 


able[FindRoot (Zeta(s] = 0, 
{s, 1/2+N[2+7+Exp(1] « (n- 11/8) /Exp(1] /LambertW{(n- 11/8) /Exp(1]), 
125} #1} , {n, 50}] 


{5-205+ 4.1347} |11} s->05+529703} | 21} {s->05 + 7.33741} | 31] {s->05 +103.726)} atl (s->05+ 242571 
S->05+ 2102001} {12} {s>05+564462 } | 2) {s->05+ 8.91041} | 32] {s->05+105-47 lal s->05+ 07507) 


wh 
f i 2 3 
{5->05+2501081} [13] s->05 45934709 | 23 {5-05 + 8473551} [33] 05 +107694 
{5->05+204081} |14/f¢>05+608389 [24l s->05-+a742s3i} [34lf6->05 1128758 [a f¢-r05 +1108) 

} | } 

} | } 

i | 

¥ 


ef 5 
{8-205 +37.58@2 I} |16]{s->05+67.0788)} | 26] {s->05 +$2.4919 |} | 36] {s->05+11420)} | 46) {s->05+ 134.757} 
{s->05+40S187|} 17} {s->05+69.5464)} | 27] s->05+94.6513)} | 37] (s>05+11627)} | 47) (s->05+ 138116 |} 
i f 


1 


{s->05+43.3271I} 18) {s->05+720672} | 28} {s->05+%5.8706 |} | 38) {s->05+1181)} |} {s->05+ 139731} 
{s->05+ 411241} 
{s->05+ 143.112 |} 


E 
4 
7] 

E 

E 

fi 


{s->05+48.0052 I} |19) {s->05+75.7047 | 2} {s->05+88312 I} | 39] {s->05+12L370)} 
{5->05+ 87781} |20]{s->05+77.148)} |) {s->05+ 103181} | 40) {s>05+12297 


Table 6: The firsts 50 complex Roots(with W-function) 
A very fast method that ends in the root very quickly and very close to the roots of ¢(s) = 0 as shown. 
5.6 Directly (from Explicit form) with the solution of this equation [8] 


We know that the equation 


tr bi . l 1 . 11 
lor ( ) + i arg (5+ 8-+it,) =n 5. 


Qn Qnre 60+ 7 8 
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Using the initial value the relation 


the explicit formula.. 


As already discussed, the function arg(¢(1/2 + 6+7-y)) oscillates around zero. At a zero it can be well- 
defined by the limit, which is generally not zero. For example, for the first Riemann zero y; = 14.1347 the 
limit 6 > Ot has value as.. 


lim arg C(1/2 +6 + é- y) = 0.157873. 


The arg ¢ term plays an important role and indeed improves the estimate of the n*” zero. We can calculate 
by Newton’s method, and we locate the first 30 imaginary part roots, of the equation ¢(1/2 +7-y) = 0, on 
the bottom of the (Table.7), and where « = Im(1/2 + yi). The process we use here does not cease to be 
approximate, although it gives us very good results quickly. The real solution to the problem, that is, the 
correct solution of the before transcendental equation, is dealt with most effectively in Chapter 5.7. 


5.7 From Explicit form with the solution of method of Periodic Radicals.. 


| Table[ 


FindRoot(x/ (2) *Log[x/ (2 m*e)] +1/x*Arg[Zeta(1/2+x«I]] =n-11/8, 
{x, 1/24+N[2s7*Exp(1] * (n-11/8) /Exp[1) / LambertW/[(n-11/8) /e], 20)}], 


{n, 1, 50}) 


1 } } {x -> 124.257} 
2| {x -> 21.0220} | 12] {x -> 56.4462} | 22| {x -> 82.9104} {x -> 127.517} 
3] {x -> 25.0109} | 13] {x -> 59.3470} | 23) {x -> 84.7355} -> 107.169} | 43] {x -> 129.579} 
Al {x -> 30.4249} | 14] {x -> 60.8318} | 24] {x -> 87.4253} -> 111.030} | 44] {x -> 131.088} 
5] {x -> 32.9351} | 15] {x -> 65.1125} | 25| {x -> 88.8091} -> 111.875} | 45] {x -> 133.498} 
6| {x -> 37.5862} | 16] {x -> 67.0798} | 26] {x -> 92.4919} -> 114.320} | 46] {x -> 134.757} 
7| {x -> 40.9187} | 17] {x -> 69.5464} | 27| {x -> 94.6513} -> 116.227} | 47] {x -> 138.116} 
8| {x ->43.3271} | 18] {x -> 72.0672} | 28] {x -> 95.8706} -> 118.791} | 48] {x -> 139.736} 
9| {x -> 48.0052} | 19] {x -> 75.7048} | 29] {x -> 98.8312} -> 121.370} | 49] {x -> 141.124} 
O| {x ->.49.7738} | 20| {x -> 77.1449} | 30} {x -> 101.318} | 40] {x -> 122.947} | 50] {x -> 143.112} 


Table 7: The firsts 50 roots of ZetaZero.(Arg() Method) 


Using the relation Explicit form.. 


tn 
20 


= 


1 
lim —a 
pig es 


tn \ 
27e 


(3 


We divide all the factors of the equation by 1/e and we take the more specialized form, according to the 
replace « = t,,/(27e) therefore we have the form 


x: log(a) + 1/(7-e)-arg(z(1/2+1-(2-7-e-2x)) =1/e- (n— 11/8) 


There are 2 ways to solve this equation or find the inverse of the expression x - log(a) or of the expression 


arg(z(1/2+1-(2-a-e-2))). Respectively the inversions are for each performance 


and 


fo(u) = Zeta* (arg~'(u)) E 


ft (u) = eFroductlog(u) 


The preferred procedure for finding a solution is the first because it is 


simpler, because the second it’s more tricky program. Following is a program in mathematica. 
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ClearAll; croots=25; 

Off[FindRoot::jsing,FindRoot::stop,FindRoot::\stol|; 

fi[x_]:=-1/(1*e)*Arg|Zeta[1/2+1*(2*2*e*x)]]+1/e*(n-11/8); 

f[x_]:=x* Log[x]+1/(a*e)*Arg|Zeta|1/2+1*(2*2*e*x)]]-1/e*(n-11/8); 

For|n=1,n<=croots,n++, 

flu_ ]:= = eee. 

xr=Nest[f[fl [#]|&,0,5]; 

FQ2=N[xr,25]; 

Print["x(",n,")="",s=N[y]/.FindRoot[h[y|== 0 ,{y,N[FQ2]},WorkingPrecision->40]; 
sl=Re[N[y1]]/.FindRoot|Zeta[1/2+I*y1]== 0, {y1,N[(2*2*E*s)]},WorkingPrecision->45], ","",.Im[N[ZetaZero[n],10]]]] 


This program is based on the reversal of x log x term, the reverse of which is associated with W - Function, 
a process that gives us a multiplicity to the Improved programming, which we take advantage in order to 
find a plurality of roots. Those with the creation of this small Program we find it quickly and simply, any 
set of complex roots we need. As we observe with this process we achieve a better approach and the roots 


x(4)=30.42487612585951257432106 | x(29)=98.83119421819369223332442 


x(13)=59.34704400260235246378215 | x(38)=118.7907828659762173229791 
x(14)=60.83177852460980984425990 | x(39)=121.3701250024206459189455 
xO2 65.11254404808160295506087 | x ele 122.9468292935525857178797 

x(20)=77.14484006887480537268266 | x(45)=133.4977372029975864501305 
x(24)=87.42527461312522666681858 | x(49)=141.1237074040211237619404 
x(25)=88.80911120763446936192582 | x(50)=143.1118458076206269771283 


Table 8: The firsts 50 roots of ZetaZero.(alog(a) Method) 


are not removed from their correct values Perhaps it is the only case of finding complex roots in a large 
area and number that we could achieve computationally. Unlike case # 5.6 which is widely used based on 
the approximate type, 


from explicit formula we may have discrepancies. That is, what is clearly seen in the calculation of complex 
roots of ¢(z) = 0 is the fact that they must treat the equation as a whole and not in parts. Finally, as 
we saw in this chapter, we fully tackled the transcendental equation as a whole and formulated without 
gaps-free algorithm with a simple process that works perfectly and efficiently. 
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6. 100% of the zeros of ¢(z) = 0 are on the critical line 
The elementary in our proof will be broken down into 5 parts. 
6.1.1 Using direct method of inverse of ¢() on functional equations. 


From the proofs of Theorems 1,2,3 we find it necessary to find the inverse of the function ¢(z) = 
¢(z) — ¢(1 — z) by the lagrange method. First of all we have to mention why we got this equation now 
and we want to examine its inverse. we saw in Theorem 3 since we proved from Th.1 that it is 1 — 1 and 
therefore the function ¢(z) is also inverted. We saw that the condition for the functional equations to hold 
is necessary in the form ...,¢é(z) = 0. But from the fact that ¢(z) is 1 — 1 since the partial functions that 
make it up are also 1—1. So ¢(z) = ¢(1—z) => z=1-—z = Re(z) =1/2. Also a very simple example is the 
fact that the series 1—* + 27-7 +37-7+...+n77 = 1771427714 37-14...+n?7! it is very simple to prove 
that take as solution the unique real z = 1/2 comparing the terms one by one. Therefore, according to 
Riemann’s functional equations, we have one critical critical line on which the non-significant 
zeros will lie. 


6.1.2 Solution of the functional Equation ¢(s), s Complex number 


For real values with Re(s) > 0 and with form 


l-s 


nq 2T (=) ¢(s) = nS TL 


2 
galt) 
¢(s) = ri) 


2 


)ca-95 


¢(1— 8) 


if we take the logarithm but or z (we are making a change in the variable s), of the two parts of equation 
[#4. Eq-set(IIT)] then we have 


cle) 6(1 = 2) = a ta@mytr (75) pr (3) => 


log[¢(z)/¢(1 — z)] = (; i :) log[m] + log Ir (- 5 *) iv (5)| | Qkeri 


The resolution analysis will be performed by the Lagrange inversion theorem, which distinguished the best 
approach for transcendental equations. Using the correlation theory and after relations with 1 from 3 groups 
fields and therefore for our case we will get: 


pi(z) = s which means that f(s) = p;'(s) = z, but with an initial value 


e atk 
a log(z) 


and total form from theory Lagrange for the root is: 


Program in mathematica 
{ Needs NumericalCalculus” | 


k := 1;q := 25; 
z= N[sin|+ 


5S ((-1/log[n])” / Gammalw-+ 1]) # 


wel 
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ND{(s’) * log[Gamma|1 — s]/2]/ Gamma[s/2]] — log[Zeta[s]/ Zeta[1 — s]])”, {s, w — 1}, 
t, WorkingPrecision + 250, Scale + .00001] 


FQ = N[z/.Sin — t, 25]; N [tFQ — 1/2* Gamma [(1—FQ)/2]/ Gamma [FQ/2] - Zeta [FQ]/ Zeta [1 — FQ]] \ 


The results for k = 1 and q = 25 limit indicator of series is 0.5 + 5.48879293259.. and with approxi- 

mation 2.22044 « 10715 + 2.664535 « 107 !°I of equation ¢. Note : The program is in language mathematica 
w-1 

and where ND is -}.-+ and s’ (is first derivative of s) = 1, also k € N°. 


Because of its limit summation tens to zero for a large n, we can roots of the functional equation to 


write it simply as: 
1 2-k-r-i 
z>n-4 : Ls EZ 
2 log(z) 


We see, therefore, that the roots of the functional equation of the ¢, are on the critical line, with Re(s) = 1/2, 
a very serious result for its complex roots. Below we see a list with k = 0,1, 2,..10 first values which are the 
roots with a very good approach of the functional equation of the ¢() and if z a zeros we have solutions: 


{z == 0.5,z == 0.5 + 5.4887929325942295i,z == 0.5 + 10.977585865188459i,z = 
= 0.5 + 16.46637879778269i,z == 0.5 + 21.955171730376918i, z = 
= 0.5 + 27.443964662971148i,z == 0.5 + 32.93275759556538i, z = 
= 0.5 + 38.42155052815961i,z == 0.5 + 43.9103434607 53836i, z = 
= 0.5 + 49.399136393348066i,z == 0.5 + 54.887929325942295:} 


If we want to see the roots of ¢(z) = 0 ie analytically for the first 12 values: 


{z == 0.5 + 14.134725141734695i,z == 0.5 + 21.022039638771556i, z = 
= 0.5 + 25.01085758014569i,z == 0.5 + 30.424876125859512i,z = 
= 0.5 + 32.93506158773919i,z == 0.5 + 37.586178158825675i, z = 
= 0.5 + 40,9187190121475i,z == 0.5 + 43.327073280915i,z = 
= 0.5 + 48.00515088116716i,z == 0.5 + 49.7738324776723i,z = 
= 0.5 + 52.970321477714464i,z == 0.5 + 56.446247697 063391} 


we can compare them to the previous values that we found of the functional equation of the ¢. Below 
is given the graphical representation of ¢- function in C for x € (0,1) and y € (—31, 31), 
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Figure 7: Plot of ¢(z) function 


we notice in the diagram Fig.7 that for the blue lines apply. Im Gee) = 0 and for the Red lines 


Zeta(1/2—I-y) 
Zeta(1/2+I-y) \ __ 
Re (aes) = 0. 


6.1.3 Solutions of the Equation ¢(s)/¢(1 — s) = 1,s is Complex number. 


A functional equation that hides all secrets of the ¢(s) = 0. Based on the analytical extension of the 
zeta function we conclude that ¢(1— s) is analytical at C — {0} and zero will be a simple pole. So we will 
have: 

lim ¢(1 — z) = co 

z—0 


lim tee) (27)* = 00 


z30 isin Fz 


ie equality applies to zo, where 
A(z) = lim, A(z) = —2isin(mz)T'(z)¢(z). 


As for the regard 4 categories roots of the equation F(s) = ¢(s) — ¢(1 —s) = 0 and s = x + yot we have a 
brief summary of the form: 


x0 = 3; Yo; if (yo = 0) 
Xo, £yoi : 

if (m0 # Frye #0) 
1- Xo, £Yol 

— 

XO, 7 

if (xo # 3,yo = 0) 
1— v0, 


Xo =1l—x = 5, if (yo =0) 


The equation F(s) is related to the Hypothesis Riemann and helps in the final proof of Theorem 3(pages 
8-10) seen before. It is a general form that also contains the non-trivial roots of ¢(s) = 0 and ¢(1 — s) =0. 
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Below we see a list with 15 first values, which are the roots with a very good approach of the equa- 
tion of the F(s) = 0 with method Newton in language Mathematica. 


Reduce[Zeta|z]==Zeta[1-z]&& Abs|z]<15,z] 


z = {—12.16278 — 2.58046 I}||z = {—12.1627878 + 2.5804649 I}||z = {—7.9909145 — 4.5105941 T}||z = 
{—7.99091453+4.51059414 I}||z = {0.5—14.13472514 I}||z = {0.5—9.666908 I}||z = {0.5—3.4362182 I}||z = 
{1/2 D||z = {0.5 + 3.4362182 I}||z = {0.5 + 9.66690805 I}||z = {0.5 + 14.1347251 T}||z = {8.990914533 — 
4.510594140 I}||z = {8.990914533 + 4.51059414069 I}||z = {13.1627870 — 2.580464 I}||z = {13.16278786 + 
2.58046497 I} 


Further, is given the graphical representation of F'(s) = ¢(s) — ¢(1 — s) functional equation. 


| a a a a a a a a a Bel a oe ae a 


-20- | 


Figure 8: Plot of ¢(z) — ¢(1 — z) function 


6.2. With general inverse random function theory. 


According to the theory of inversion of functions that help us to solve them as equations, it completely 
guides us to a single solution i.e. Re(z) = 1/2 for ¢(z) = 0. In this article, we seek to develop a formula 
for an inverse Riemann zeta function w = ¢(s) we have an inverse function s = ¢~!(w) and whics implies 
that ¢~!(¢(s)) = s A ¢-+(¢(w)) = w for real and complex Domain w and s. The presented method can 
also recursively compute these multiple solutions on other branches, but as we will find, the computational 
requirements become very high and start exceeding the limitations of the test computer, so we will primarily 
focus on the principal solution. We develop a recursive formula for an inverse Riemann zeta function as: 


£4)" 


s—0 k=1 5k 


5 =6Mw) = lim m |— gy ae og [(6(8) - w) (6 — 0) 


m—0oo (2m = 


Where S,, is the multi-value for which w = ¢(s). 
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6.2.1. Theorem 4. To prove that for the formula [26] 


1 


ae log |f(z)|z+0 = Z(m) — P(m) 


that valid for a positive integer variable m > 1. This formula generates Z(m) — P(m) over all zeros and 
poles in the whole complex plane, as opposed to being enclosed in some contour. 


Proof. If we model an analytic function f(z) having simple zeros and poles by admitting a factorization of 
the form 


Np 4 
(=r (-=)-T (cr 
n=1 n=1 
Where g(z) is a component not any zeros or poles, then 
Np 
log|f(z)] = log[g(z)] + s log (1 = =) ~ Slog (1 = =) 


From theorem Residue if i take 


And with use Taylor Theorem and |z| < 1 we take in final 


[oe) [o.e) gk 
log[f(z)] = logg(z)] — S> Z(k ae S> P(k) = 


=1 k=1 


x 


From this form, we can now extract Z(m) — P(m) by the nth order differentiation as 


Capa 0 = Aon) ~PCm) 


6.2.2. Theorem 5. If Z,, is a set of positive real numbers ordered such that 0 < z1 < z2 <...< Z,, and 
so on, then the recurrence relation for the nth+1 term is 


‘i 1 —1/s 
Zn41 = lim (209 - z| 
k=1 


Proof. 


The generalized zeta series over zeros Z, of a function 


z 


N 
1 1 1 1 
Z@)= 05 = tate 
n 1 


s Ss 
Z5 Z3 


n=1 


and let us also assume that the zeros are positive, real, and ordered from smallest to largest such that 
O0< 2 <Z2<...< 2%, then the asymptotic relation holds 


1 


. =O (2) > Ole 8s 0 
n+1 
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Final we take 


1 —1/s 
—1/s _ =. (as 
aa = fim 2(6))-N, 2a = Jim (26) 5) so 
i 1 —1/s 
to = am, (20)-S 2) 
k=1 “k 


thus all zeros up to the nth order must be known in order to generate the nth+1 zero. Exist the Voros’s 
closed-form formula for depends on (RH), and the formula which we can express the zeta and beta terms 
in terms of a Hurwitz zeta function, and then substituting the the Voros’s closed-form formula we obtain 
another formula for non-trivial zeros a for in terms of the Zeta .If the limit converges to 1/2, it would imply 
(RH) is correct. The formula of Voros’s is especially for real part (of pune the non-trivial zeros) of 


¢(s) = 0 is 


1 1 
R (Pint) =o1 = lim (5226) a 5Zm(2m)) + 


m—>co 


(= atten (5) am (2m3))) | =3 


m R(p1, nt) 

15 | 0.4730925331369 
20 | 0.4898729067547 
25 | 0.4993065936936 
50. | 0.5000000028549 
100 | 0.4999999999999 
150 | 0.5000000000000 | Table 9: Table Voro’s. Approximation of 
200 | 0.4999999999999 Re(z) of z() function 


V.roots 


The computation of the real part of the first non-trivial zero for m = 200 places (13 digits) using this 
relationship Vroots and with 200 iterations we achieve a very satisfactory result. It is of course obvious that 
it approximates by the direct inverse of ¢() , the real part of the non-trivial roots to 1/2. 


6.3. Statistical evidence 


Here, we formulate two of the functional equations from E-q.Set (Theorem 2, page 7) 


(1 — s)/¢(s) 
¢(s)/¢(1 — 8) 
We look at each one equation individually in order to identify the set of values that we want each time. 


For the first equation and for real values with Re(s) > 0 and by taking the logarithm of two sides of the 
equation [5], we have... 


2(2m)~* Cos(ms/2)T(s), s > 0 


2(27)*—' Sin(rs/2)T(1—s), Re<1 


log|2| 


real s we will have | Re(s) > = 0.3771]. This is the lower bound, which gives us the first ¢(s) 


of Riemann’s Zeta Function. 


For the second equation, for real values if Re(s) < 1 and by taking the logarithm of the two parts of 
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the equation, we will have: 


log|z 


In the following, because we need real s we will take | Re(s) < = 0.6228 


average value of the upper lower bound is Re(s) = 1/2 because from that is apply 


_log[2] + log{z] _ 


Re(s) = 1/2: TP = 1/2 


6.4. The critical line also follows from the definition of A(N) 


Will define the function \(V) as follows: 


A(N) =14+2754+384...4N5 


The next step is to perform a mathematical analysis on the functions of after allocating N value by 
equalizing it with P,i1, because the values of X function in mainly depends on the prime numbers. In 
this way, the equations that represent the imaginary part of zeros which we are going to have will be as a 
function of the prime numbers. But the equations include all the natural numbers, so we’ll generalize the 
solutions which depend on the prime numbers after obtaining all natural numbers with the same values 
which approach to those prime numbers. 


Now we get the equation... 


1—s)p,* 1 1 1 1 
258-6) (-(-ah) 
Pn4+i 2 3 Pn-1/ Pn 
With replacing 1 — s by s and vice versa in the same equation, we get the 
— s)ps-} 1 1 1 1 
ee are 
Pr+1 2 3 Pn-1/ Pn 


By dividing equations the previous two, we get rid of the two factors s and (1 — s) which cause the 
emergence of the classical zeros of eta function 


Ss —s 2s—1 
reat ae Te al 
Pn Pn4i Pn 


s—1/2 — ago 

py Pn+1 ° Pn+1 (s-4) a 
=F s—1/2 =0 

Pn+i Pn Pn 
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And therefore 


2% 


=> an fe 1/2) In (22) =0 


Pr 
=> -i(s- 1/2) in ( | = +kr,k =1,2,3,... 


ie" 1/2yin(Fptt) i(s—1/2) n( 2p) : = 


ee ey a ee 


2 In (25#*) : 


(6.4.1) 


this last equation (6.4.I) proves the truth of the Hypothesis Riemann. 


6.4.1. A very important inequality arises from the 3rd functional equation for Number Theory 


From the 3rd Solve of the functional equations equation of ¢(s)(page 10, III) we have for real values 
with Re(s) > 0 where s is Complex number and with form 


x 


a-iL (5) ¢(s) =a ST (- 5 *) ¢(l—s)=> 


(oy =a De _ 5) 


if we take the logarithm but or z (we are making a change in the variable s), of the two parts of equation 


then we have 
=) (5) 


lou(e(2)/€(1~ 2)) = (5 ~ 2) loe() +108 (r (+ 


¢(2)/¢(1 — 2) = 77742(2m)? OP ( 


*) iv (=) + 2kri) 


The resolution analysis will be performed by the Lagrange inversion theorem, which distinguished the best 
approach for transcendental equations. Using the correlation theory and after relations with 1 from 3 groups 
fields and therefore for our case we will get 


pi(z) =s which means that f(s) = p;'(s) =z but with an initial value 
1 2-k-m-i 
Sn 


2 log(z) 


And total from Lagrange inversion theorem for the root Z is if q > co 
2 = Sint « —1/log|m])”/Gamma|[w + 1}) 


(2s (log|(Gamma|(1 — s)/2]) * (Gamma|[z/2])] — log[(Zeta|a]) (Zeta[1 — x])])”’ 


Because of its limit summation tens to zero for a large n, we can roots of the functional equation to write 
it simply as: 


1 2-k-a-i 
+ 


> Z 
2 log() ? 
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We see, therefore, that the roots of the functional equation of the ¢, are on the critical line, with Re(s) = 1/2, 
a very serious result for its complex roots. Below we see a list with k = 0,1, 2,..10 first values, which are 
the roots with a very good approach of the functional equation of the ¢ and if z we have solutions: 


{z == 0.5, z = = 0.5 + 5.4887929325942295t, z == 0.5 + 10.9775858651884591, z = 
= 0.5 + 16.466378797782697, 2 == 0.5 + 21.955171730376918:7, z = 
= 0.5 + 27.4439646629711482, z == 0.5 + 32.93275759556538:, z = 
= 0.5 + 38.421550528159612, z == 0.5 + 43.9103434607538362, z = 
= 0.5 + 49.3991363933480662, z == 0.5 + 54.887929325942295t} 


we can compare them to the previous values that we found of the functional equation of the ¢(z) equation. 
But from the second proof the non-trivial zeros z has a lower value than the first one. This means 
that that is: z < s(from 6.4.1) that is 


2KTt KT 


< 
log(7) log (4 


n 1 n 
=> log (2 +) < =log(r) > 1 < Va 
) Pn 2 Pn 


Sure it is obvious 


Bee ts ie OEE a (6.4.2.1) 
Pn Pn 


The ratio of the largest to the smallest consecutive prime has an upper bound of the real number ./7. A 
relation that enters into many prime number conjectures. 
6.5. From the inverse of the function ¢() 


First we need to formulate the Riemann functional equation. It is necessary to define the function €(s) 
1 —s/2 1 —s/2 I 
&(s) := 55s -—)arerrr 3° ¢(s) =(s—-1)n */°T 5s +1) ) ¢(s) 
Where we have used the identity ssl ($s) =. ($s + 1). 
We use the following theorem, 


Theorem 6. The function € has an analytic continuation to C. For continuation we have €(1 — s) = &(s) 
for s EC. 


6.5.2. Corollary 2. The function ¢ an analytic continuation to C\{1} with a simple pole with Residue 1 
ats=1. 


For this continuation we have 


¢(1 — 8) = 2'-*x~* cos (5"s) I'(s)¢(s) for s € C\{0, 1} 


Proof. We define the analytic continuation of ¢ by 


a) 78/2 1 
&(s)n°" Sr sti) 
s—l1 


¢(s) = 
1/T is analytic on C and the other function in the numerator are the analytic on C.Hence ¢ is analytic on 
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C\{1}. The analytical continuation stated here coincides with that stated by the well-known theory {s in 
C:Re s > 0}\{1} since analytic continuations to connected sets are uniquely determined. Hence ¢(s) has 
a simple pole with residue 1 at s = 1. We derive the functional equation. By Theorem 6.5.1 we have, for s 


in C\{1}, 


_ (1s) _ &(s) 
$(1—s)(—s)m-G-9)/2P ($(1—s))  $8(s —1)m-G-9)2T (F(1 — 8) 
1, eo —s/2 1 
ee ee ay OPO 


58(s —1)7 —(l- 20 (4 (1- 


Now we have 


This implies this Corollary. 


6.5.3. Direct solution of the equation ¢(q * z) =0,qeER 


Keeping the basic logic of the corollary 6.5.2., declare the definition s + q-z,q € Z who is analytic 
on C, for every fixed q € R,q #0. Then is obvious. According Corollary 6.5.2. we have 


1 
¢(1—q-z) =2'- 97 q7—4? cos (5" -q: 2) T'(q:z)¢(q-z) for q-z € C\{0, 1}. 
But according to theorems 1 & 3 the function ¢() is 1-1 in the critical strip, which means that: 
¢(1—q-z) = ¢(q-z) and therefore 


1 
z= =,q40,q ER (a) 
2q 


This relationship always determines the position of the critical line, and of course atom is directly related 
to the ¢() function as we proved in Theorem 3. Let us now define a function ZetaZero() which gives us the 
zeroes of ¢(s) = 0 as defined in sections {5.4, 5.5, 5.6, 5.7}. 


We therefore have the relation s = ZetaZero(£k),k € N,k > 1. Therefore 


s = Re(ZetaZero(k)) + Im(ZetaZero)(k) -i, i is imaginary unit, k € NT (b) 


As we know that for the simple function ¢(s) = 0 the 


1 
Re(s) = Re(ZetaZero(k)) = ak e€ Nt 
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From relations (a,b) we have the relation for non trivial zeroes, 


1 1 
Zneroes = — * Re(ZetaZero(k)) + — - Im(Zeta Zero(k))-ie 
q q 


1 1 
Zreroes = =~ £ — - Im(ZetaZero (k)) «i, 
5g tS Im( 0s) 
iis imaginary unit,k € N*,q € R, where q 4 0. (c) 


from the last relation (c) it follows that the Critical line will be given by the relation 


- Re(ZetaZero(k)) © 


1 
Zeline = 
q 
1 
Zeline = ma € Nt, qeER, where q x 0. (d) 


Finaly for the case of trivial zeroes we have the relation, 


Zzeroes = Jk Nt, q R, where q x 0. (e) 


With this analysis we have essentially achieved the general solution of the transcendental equation ¢(q:z) = 
0,q € R, and we observe that for positive values of q, the critical lines, infinite in number, run through the 
interval (0,1) for three cases (1&2&3), 


1 1 1 
1. 0 < Zeline < = then 0 < — < = therefore q > 1 
2 2q 2 


1 1 1 1 
2. 5 < Zeline < 1 then 5 < 2a < 1 therefore 3 <q<l 


1 1 
3. Zeline > 1 then — > 1 therefore 0 <q < = 
2q 2 


Similarly for the interval (—1,0) we will have the relations (4&5&6), 


1 1 1 
4. ~ < Zeline < 0 then — = < — <0 therefore g < —1 
2 2) 2G 


1 1 1 1 
5. —1<Zoline < 5 then —1< aq < 5 therefore -l1<q< 5 


1 1 
6. Zetine < —1 then — < —1 therefore —-=<q<0 
2q 2 


6.5.4. Two Characteristic examples 
For the Riemann hypothesis the most important thing we are interested in is the interval (0,1). The 
transcendental equation z(q* z) = 0 as we have shown has infinite critical lines in this interval, which 


depend exclusively on ato q. 


i) I take for begin the mean value of interval (0, 1/2) and we have: 
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1 1 
1. Zc. line =4=3309=2 


non trivials zeroes 

Ze oee al ZetaZero (+m)),m = 1+5 

Zzeroes = {0.25 + 7.067367, 0.25 + 10.5117, 0.25 + 12.50547, 
0.25 + 15.21247, 0.25 + 16.46757} 


trivials zeroes 
—_ m e . 
3. Babeses = 2x Als 


Zzeroes =f Le 2s, 3., 4., 5.} 


ii) Also the mean value of interval (1/2, 1) is 3/4. Therefore we get the relations for this case: 


1. Zc. line =f=%049=3 

non trivials zeroes 

Paes Gl ak ZetaZero (+m)),m=1+5 

Zgeroes = {0.75 + 21.2021T, 0.75 + 31.5331, 0.75 + 37.51631 
0.75 + 45.6373I, 0.75 + 49.4016I} 


trivials zeroes 
3. Zzeroes = —2 X = ,m =1+5 
Zzeroes = { 3., 6., 9., 12., 15.} 


6.6. Deniers of the Riemann Hypothesis. 


The most characteristic case that can be considered in relation to ¢(z) = 0 is the Davenport-Heilbronn 
function (introduced by Titchmarsh), is a linear combination of the two L -functions with a complex char- 
acter mod 5, with a functional equation of L -function type but for which the analogue of the Riemann 
hypothesis fails. In this lecture, we study the Moebius inversion for functions of this type and show how its 
behavior is related to the distribution of zeros in the half-plane of absolute convergence. In general, however, 
these cases do not refer to ¢(z) = 0 but to the L -function. The most important element in the negators 
of the Riemann hypothesis to consider is whether it is consistent with the Riemann equation ¢(z) = 0. 
Therefore, we will not consider cases that do not fall within our hypothesis, as defined by Riemann. 


6.6.1. Introduction of the denial conjecture through control one case of ¢(6, z) = 0. 


However, there is one interesting case announced by the Enrico Bombieri .In a video by Enrico Bombieri[23] 
it is clear that there is a root that belongs to the 


€(s, (m2? +5-n2)) =0 


The root is s = 0.93296997 + 15.668249531*I and is the prime for m = n = 1. Thus the equation becomes 
¢(s,6) = 0 as he mentions. We have to prove that it is not the root of the equation. 


True that ¢(s,6) = —0.101884 — 0.005983947 approximately. In contrast, the ¢(6,z) case gives a value 
to zero. 
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In more detail approximately, 


C(6,z) = 2.93166 - 1078 + 2.10446 - 107" - I 


But if we do a simulation with the equation ¢ (q*z) = 0 then we need g = 1/(2(Re(z))) = 1/ (2*0.93296997) 
i.e q = 0.535923 so we will have another root z = 0.93297 + 26.3745I approximately, for which 


C(q: z) = —2.46975 - 107 1° + 5.81423 107"° - I 


a very good approximation and correct this time. To verify the fact that this root is a false root we will use 2 
methods. One is with the Newton method locally and the other is to solve also locally the transcendental 
equation ¢(6, z) i.e. near the real root of the complex Re(s) = 0.93296997. If in this zone we do not find 
such a value this root will no longer exist. 


6.6.2. Approximate proof of the Newton method for checking ¢(6,z). 


So assuming that the root is actually sl = 0.93296997 — 15.668249531*I, we set it as the initial root in 
FindRoot|Zeta [6, z], {z,s1}, WorkingPrecision — 50] with a dynamic iteration 50 times. We imme- 
diately see that it locally diverges to infinity and gives us a value 


zZ — 3.5898 - 10" — 1.29617 - 10°I 


This means that it failed near there to approach any root that would zero the equation ¢(6, z). In contrast, 
by the same method we observe that the local root s1 = 0.93297 + 26.3745I for the equation ¢ (q*z) very 
quickly locates root i.e. using the command FindRoot|Zeta [q*z] , {z,s1}, WorkingPrecision — > 50] 
with q = 0.535923 then z = 0.93296997 + 26.3745478187 with approximation 


¢(q-z) © 2.46975 - 107° — 1.29617- 107”. 


Therefore the root s1=0.93296997-I*15.668249531, failed the first test to confirm that it is a root. 


6.6.3. Direct proof of ¢(6,z) = 0. 


In mathematics, the polygamma function of order m is a meromorphic function on the complex 
numbers C, defined as the (m+ 1) th derivative of the logarithm of the gamma function. 


It satisfies the Reflection relation: 


(—1)"bGm, 1 — 2) — v(m, z) = m5 cot(m - z) (6.6.3.1) 
The polygamma function has the series representation 
i Backs 
v(m, z) = ml) GP et 5 aa (6.6.3.2) 


which holds for integer values of m > 0 and any complex z not equal to a negative integer. This represen- 
tation can be written more compactly in terms of the Hurwitz zeta function as 


(m,z) = (—1)™** m!¢(m + 1,z) (6.6.3.3) 
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Alternately, the Hurwitz zeta can be understood to generalize the polygamma to arbitrary, non-integer 
order. If we generally define ¢(m+1,z) = f(z) we can with improved programming solve this equation based 
on the 2 previous equations (6.6.3.1) & (6.6.3.3). In this case we ask to find the roots when ¢(m+1,z) =0 
therefore f(z) = 0 with method Periodic Radicals. 


In this case we are interested in the fact that ¢(6,z) = 0 and therefore ~(5,z) = 0 from relation (6.6.3.3). 
From (6.6.3.1) we have 


5 
(CDG iS s06 esa _ eee) (6.6.3.4) 
But because 0(5,z) = 0 the resulting final relation will be 
ap 
—(5,1—z) =7- —scot(m- z) (6.6.3.5) 
0z° 


Performing the process of finding roots according to the Generalized Root Existence Theorem (of an equa- 
tion) [24] we find the complementary function after solving the form 


w(5,1—z)+n°- (-16- (cot(az))* - (esc(mz))? — 88 - (cot(mz))? - (cse(wz))* — 16 - (csc(rz))°) =0 
(6.6.3.6) 


in final form 


(5,1 —z)/n®° + (-16- (cot(mz))? - (ese(mz))” — 88 - (cot(mz))? - (ese(mz))* — 16 - (ese(mz))*) = 0 
(6.6.3.7) 


To solve such a transcendental equation, we will use Periodic Radicals, which means we will split the equa- 
tion into 2 parts. We solve the 2 nd part because it’s a trigonometric equation and we have, 


(—16- (cot(mz))* - (ese(mz))? — 88 - (cot(z))? - (esc(mz))* — 16- (cse(mz))°) =0 


whose has in general 12 solutions j := +1;m := 1 + 6; as below, 


=] 


(i +2 ArcTan | VRoot [15 + 30#1 + 17H 12 + (44 8u) #13 + ITH! + 20K + 15H 18k, m| ae ark) 


23S 


A radical relation which gives us the inverse of the 2nd part of equation 6.6.3.7.The whole secret is based 
on the parameter k of the trigonometric equation, which in turn will give us multiple different roots, which 
solve our basic equation ¢(6,z) = 0. The total cases we consider are 2X6=12, and of these we keep the 
values where 0 < Re(z) < 1. With these zeros we make the following table of the program below. 
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Then the matrix program in mathematica witn infinity Periodic Radicals is, 


ClearAll; j:= +4;m:=1+6;n:=5;d:=5; Off[FindRoot::Istol,FindRoot::precw,Divide:nfy, 


Infinity::indet,General::cvmit,ReplaceAll::reps,FindRoot::jsing]; 

h2[x_]:=-(0+(-1)*(n+1)*PolyGammaf[n, 1-(x)])*1/ (26); 

h[x_]:=0+(-1)\(n+1)*PolyGammaf[n, 1-x]+a° (-16 Cot[z x]* Csc[a x]?-88 Cot[z x]? Csc[a x]*-16 Csc[m x]9); 
h3[x_]:=N[PolyGamma|[n,x],50]; 

h4[x_]:=N[Zeta[6,x],50]; 

For[k=-d,k<= 0,k++, 


flu ]=I/n (j*2 AreTanfV Root [15 + 3041+ 1741 + (4+ 8u) #1° + 17417 + 3941 + 154° &, m] }22h);: 
xr=Nest[f[h2[#]]&,-1/10%7,1]; 
FQ2=N[xr,40]; 
Print["x(",k,")=",s=N[y,40]/.FindRoot[h[y]== 0, {y,FQ2}, 
WorkingPrecision->25],",",N[h3[s]]]] 


The final table of results is given below, 


k/j=+1 m=1 m=2 m=3 m=4 m=5 m=6 
5 -9.4974772-2.4579566i | -9.4974772+2.4579566i | -9.5000000-0.500433 i} -9.5000000+0.500433 i} -9.499999-0.134109i | -9.499999-0.134109 i 
-4 -7.4923956-2.4581002i | -7.4923956+2.4581002i | -7.5000000-0.500433 i] -7.5000000+0.500433 i} -7.499999-0.134109i | -7.499999-0.134109 i 


3 
2 
A 
0 
1 
kj | med meme mes me 
5 
4 
3 
2 
A 
0 
1 


Table 9: Zeros of equation ¢(6, z) = 0 


As we observe Re(z) < 0 and k < 0, for all roots. If we solve the equation for k > 0 we obtain exponential 
roots, i.e. large roots, were simply omitted. Therefore there is no solution for Re(z) in the interval (0,1) as 
a critical line. 


6.6.4. Direct proof of Epstein zeta equations ¢(Q(m,n),s) = 0 and the contradiction of the 
deniers of R-H. 


In 1935 Potter and Titchmarsh for equation ¢(Q, z) proved has infinite zeros on the critical line and they 
also reported that a numerical calculation had shown that in the case, the Epstein function zeta has a zero 
on the critical strip but not on the critical line. This put an end to the speculation that the functional 
equation was the reason for the prevalence of a Riemann hypothesis. The deniers state ”It is widely accepted 
that the generalized Riemann hypothesis is valid, with strong numerical proofs supporting it, but the proof 
remains long elusive.” The Epstein zeta functions have the form of a double sum 
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aS r( (pi.P Guan (pi, P2) =a- pi+b- Pipa +c: pa. D= b* — 4-a-c},{p1,p2,a,b,c € Z+} 
P1,p240 i 


To complete the functions related to the ¢-function we have to deal with the ¢(4,z),¢(5, z) which the 
deniers mention. We consider the sum: 


n 


(Os) = oS 


. \2)° 
P1,p2 Pipes d*) 


where the sum over the integers p; and p2 runs over all integer pairs, apart from (0,0), as indicated by the 
superscript prime. The quantity Q corresponds to the period ratio of the rectangular lattice, and s is an 
arbitrary complex number. For Q? an integer, this is an Epstein zeta function, but for Q? non-integer we 
will refer to it as a lattice sum over the rectangular lattice. 


love) s—1/2+n 
K(n,m;s8;A) = 7” S- (oe | K5-1/24m (20p1p2d) . 


Pi,p2=1 1 


For X >= 1 and the (possibly complex) numbers small in magnitude, such sums converge rapidly, facili- 
tating numerical evaluations. (The sum gives accurate answers as soon as the argument of the MacDonald 
function exceeds the modulus of its order by a factor of 1.3 or so.) The double sums satisfy the following 
symmetry relation, obtained by interchanging pi and pg in the definition: 


K(n, —m; 8; A) = K(n,m;1— 5; A). 


Finally we get the relation for [25] 


n 


So(A, 8) =6(Q,8) = So ( : 3 Q=pi tps: 


2: 2 2 
P1,p2 Py + 2 x A”) 


200) Tis= T2)ces 1). 81° yak 
So(A, 8) = +2/r T(s)A Tisai (0.0:5:5) : 


From the general relations below we will derive the final relations of interest for the function z(Q,s) : 


So(1, 8) = 4¢(s)L_a(s), So(W2, 8) = 2¢(s)L_s(s) 
So(V3, 8) = 2 (1 — 2*-**) ¢(s)L_3(s), So(W4, 8) = 2 (1 — 27° + 2'-?*) C(s) L_a(s) 
So(V5, 8) = ¢(s)L—20(s) + L—4(s)L45(s), So(V6, s) = ¢(s)L—24(s) + L-3(s)L48(s) 
So(V7, 8) = 2(1— 21-8 + 21°) ¢(s)L_7(s). 


From the zeroing of the first part of the transcendental representations for the cases where the values for » 
are \ = /3,\ = V4, = V7 and will we take: 


So(V3,s):8 = 1 (1 “ ern) 


1 éarctan/7 = 2nzi 
So(v4,9):9= 54 —T oo + 
1 iT 2nTi 


So(v7, 8) 8 = 5 alao ae 


We note that the critical line 1/2 is technically constructed by solving the representations for these cases 
and, of course, a number of infinite imaginary roots on the critical line. And in the case of the Epstein 
Function we have points on the critical line. The only way to refute any negation for the z(Q, z) functions is 
to solve completely the 3 transcendental equations ¢(4,z) = 0,¢(5,z) = 0,¢(8,z) = 0 and check the results 
if any cases belong on the critical line. 


Solving the other cases of Epstein functional equation approximately with the Periodic Radicals method. 


If we assume that 


Q(m,n) = am? +bmn-+ en? 
= a(m+nr)(m+ nr) 
= alm+nrt/? 

= 1 


in addition, the following relationships apply 


b+ Vb? — 4ac 
—— 
2a 
r=, y=, TH=XU+Y= b+ivD 
_b 
Toe T= = 
T= 
a 
separating the term with n = 0 and we have 
Die eek OP 1 
Z(s) = Ss mes t as S- S- |m + nrp23"° > 1 (6.6.4.1) 
m=1 n=1m=—oco 


which is done 


= 1 T(s— 3) va Ar ~ fad 
S- oe = n2e—Ty2—IT (5) + nae 1y2s-TP (5) os (mnry)* 2 cos(2mnmx)K,_ 1 (2mnzy) 
m=1 


m=—Cco 


with use the Bessel function 


ie eae 1 °° cos yt 1 
Ki(y) =~ -) Piv4 | ~ dt, >0, Rw) >-=H, 
w) vn () ( >) 0 (14+¢2)"t2 ; OP 5 


and finally we get the equation for the Epstein function in the form 


oe) 


L —s 4 S65, OO 7 
) V7 e095 1)4 -_ Ts) a Py ae So (mn)*-? cos(2mnrx)K,_1(2mnry) 
n=1 


em -s. 1-268 (s- 
Z(s) =2a~5 s a7 sy—2s 2 
(5) = 20-*¢(2s) + 2amty- 


i 
2 


m=1 


(6.6.4.1) 
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It also turns out that is valid 


(7) a (2) : I(1—s)Z(1—s) (6.6.4.1I1) 


which is the functional equation of the Epstein function. But as we proved [Theorem 1, 2,3 pages 1-10] with 
the function ¢( ) that it is 1—1, the same is proved here by (6.6.4.3), which follows because Z(s) = Z(1 —s) 
ie. s=1—s=>s = 1/2 for the Real part of the variable s. But we can also prove it from equation (6.6.4.IT) 
if we consider the third part of the sum with a Bessel function as a small loss in the final result and because it 
is always imaginary, keeping the first 2 parts and obtain equation (6.6.4.IV) and work on it to find the roots. 


Z(s) = 2a-*C(28) + 2a7*y!-24~ C a VT 5-1) (6.6.4.1V) 


Thus we secure the real part which will show us the critical line and indeed uniquely. With the method 
we follow [24] (Generalized Theorem -Method Periodic Radicals) we find the generator (inverse function) 
of a part of the equation and solve the equation as a whole. For this equation a good generator is the term 
y!-28, The analysis that will result after logarithm will give us multiple different roots, which is exactly 
what we have in such a case. Below we give the relationship of the generator, 


yl 2) =u>ye=log(u)+2-k-i-7 
(6.6.4.V) 


1 log(u)+2-k-i-7 
2 (age ) 


Then the matrix program in mathematica witn infinity Periodic Radicals is, 


Clear All: 

d:=10:a=1 ;c={2,3,4,5,6,7,8,9}:y=(4 *a*c)\(1/2)/(2a); 

Off[FindRoot::Istol.FindRoot::precw,Power::infy Infinity: :indet, 
General::cvmit,ReplaceAll::reps]: 

h2[x_]:=-(2*a(-x)*Zeta[2x])/(Gamma[x-1/2]*1(1/2)/Gamma[x] *Zeta[2x-1]): 

h[x_]:=2*a’(-x)*Zeta[2x]+2*a(-x)*y\(1-2x)*Gamma|[x-1/2]*2(1/2)/Gamma[x]*Zeta[2x-1]; 

Forfk=-d.k<=d,.k++, 

flu_]:=N[1/2-(2 Ia k+Logfu])/(2 Log[y])]: 

xr=Nest[f[h2[#]]&.-1/10.5]; 

FQ2=N[xr.40]; 

s=N[y1.40]/.FindRoot[h[y1 ]== 0. {y1.FQ2}.WorkingPrecision->40]: 

Print["x(".k,")}=".s1=N[y1.40]/.FindRoot[h[y1 ]== 0. {yl.s}. WorkingPrecision->60].".".N[h[s1]]]] 


By running the program, all the roots with real part equal 1/2 are obtained. Here is a table with partial 
roots for A? = 2,3,4,5,6, 7, 8,9. 


Let’s look in detail at the actual zeros per case in Table 10. 
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L x(k) | c= 2faev2_ | ca3fdevs ca afevg | ceSfevs 
[x(-5) | _0.5+47.736836363i | 0.5+30.17305559254i| _ 0.5+24.41410967617i_| _0.5+21.050437854642 i | 
[x(-3) | _0.5+30.122983610i_| _0.5+19.4730129592i | _0.5+15.65382790863i_| _0.5+13.403898708822 i | 
[x(-1) | _0.5+7.1448574376i_| _0.5+8.17903931974i | _0.5+6.015761092906i_| _0.5+5.3257295685252 i | 
| x(0) | __0.5+/-3.603666406i_| 0.5+/-2.26693469535i| 0.5 +/-1.74489943291i | _0.5 +/-1.74489943291 i | 
| x(1) |__0.5-7.1448574376i_| 0.5-8.17903931974i | _ 0.5-6.015761092906i_| _0.5-5.3257295685252i_| 
| x(2) | _0.5-21.178130885i_| 0.5-13.6950229708i | _ 0.5-11.12944087760i_| _0.5-9.1902438527654i_| 
| x(3) |__0.5-30.122983610i | 0.5-19.4730129592i | _ 0.5-15.65382790863i_| _ 0.5-13.403898708822_| 
| x(4) |__0.5-39.264494014i_| 0.5-24.54972699805i | _0.5-19.69500841582i_| _ 0.5-17.320998857098 _| 
| x(5) |__0.5-47.736836363i | 0.5-30.17305559254i | __0.5-24.41410967617i_| _0.5-21.050437854642i_| 
x(k) | ce6/aeve | ce7fAev7_ | c=Bfaev8 | ce9fAev 
[x(-5) | _0.5+19.2325630379i | _0.5+17.4705910431i | _0.5+16.023178166 i_| _0.5+12.9585161927i_| 
[x(-3) | _0.5+11.7218046027i | _0.5+11.1453755465i | _ 0.5+9.8784324824i__| _0.5+9.50684144580i_| 
[x(-2) | _0.5+8.44537446121i | _0.5+7.93099697246i | _0.5+7.5982447947i | _0.5+6.48168939942i_| 
[x(-1) | _0.5+4.83435490426i | _0.5+4.47572646118i | _0.5+4.2020221468i | _0.5+3.98521413380i_| 
| x(0) | _0.5+/-1.2680423210i | _0.5+/-1.1330903582i | _0.5+/-1.030068120i_ | _0.5+/-0.9478863402i_| 
| x(1) | _0.5-4.83435490426i | 0.5-4.47572646118i | 0.5-4.2020221468i__| _0.5-3.98521413380i_| 
| x(2) |_0.5-8.44537446121i | _0.5-7.93099697246i | _0.5-7.5982447947i__| _0.5-6.48168939942i_| 
| x(3) |_0.5-11.7218046027i | _0.5-11.1453755465i_| _0.5-9.8784324824i | _0.5-9.50684144580i_| 
| x(4) | _0.5-15.6402085885i_| _0.5-14.0500485667 i 
| x(5) |_0.5-19.2325630379i | 0.5-17.4705910431i | _0.5-16.023178166i | _0.5-12.9585161927i_| 


Table 10: Zeros of equation Z(z) = 0 for \? = 2,3, 4,5, 6,7, 8,9 


6.6.5. Direct proof of ¢(4, z) = 0. 


The whole procedure follows from chapter 6.3 and equation 6.3.3.1 It satisfies the Reflection relation: 


re. 


(-1)"v(n,1—- z)- W(n,z) =7- aa cot(7 - z) (6.7.1) 


this equation can be written making the first part positive as 


n 


(-1)"**u(n, 1-2) + ¥(,z) = —7- _ cot(m - z) (6.7.11) 


But since we associate the polygamma function and Hurwitz zeta function arising 


(n,z) = (—1)"P*nl¢(n + 1,2) (6.7.IIl) 


As a complementary relation we take the n -th derivative of cot(7z) which will be 


yn 


—n- ~ cot(m +z) = m4 (4. cot(z)? csc(rz)” + 2 cse(mz)*) (6.7.IV) 
Z 


We come therefore to the final system 


(—1)"*!h(n, 1 — z) + ¥(n, z) = 14 (4- cot(z)? ese(mz)? + 2 csc(z)*) 
w(n, z) = (-1)"t 1 nl¢(n + 1, z) (6.7.V) 
¢(n+1,z) =0 
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Which takes the final form 


(—1)"*1 h(n, 1 — z) = 14 (4; cot(mz)? csc(rz)? + 2esc(mz)*) 
w(n,z) = (-1)"*!nl¢(n + 1, z) = 0 (6.7. V1) 
C(n + 1,2) =0 


The whole philosophy is therefore that if we admit that we zero ¢(n + 1,z) then automatically w(n, z) 
will also be zero. If we know a method that can solve 6.7.II we can also solve ¢(n + 1,z) = 0. The method 
we follow here is the Infinite Periodic Radicals, which is considered here the fastest and safest, because 
we have high polynomial powers resulting from the multiple derivation of cot(z) of the equation (6.7.II). 
In our case n = 3 to solve ¢(4,z) = 0. 


We therefore have the following program, 


ClearAll; 

d:=53n:=33m:=14;]:=2;j=-/+1; 

Off[FindRoot::jsing,FindRoot::precw,Power::infy, Infinity: :indet, 
General::cvmit,ReplaceAll::reps,Divide::infy]; 

h2[x_]:=(1/10“1*(n!)/(-1)(n+1)+(-1)(n+1)* PolyGammaf[n, l-x])/z? ; 

h[x_]:=1/10“1*(n!)/(-1)(n+1)+(-1)(n+1)*PolyGammafn, 1-x]-2* (4 Cot[z x]? Csc[ x]?+2 Csc[x x]*); 

h3[x_]:=(-1)*(n+1)*PolyGamma[n,x]/n!-1/10“; 

h4[x_]:=Zeta[n+1,x]-1/1041; 

hS5[x_]:=Zeta [n+1,x]; 

For[k=-d,k< 0,k++, 


f[u_]:=N[(j *9 ArcTan[V Root [3 +4#1+ (2-8u) #17 +41? + 3#1* &, m] }2 1 k)/n]; 


xr=Nest|f[h2[#]]&,-1/100,1]; 

FQ2=N[xr,20]; 

sl=N[y,20]/.FindRoot[h[y]== 0, {y,FQ2}, 
WorkingPrecision->30]; 

Print["x(",k,")}=",s=N[y,20]/.FindRoot[h[y]== 0, {y,s1}, 
WorkingPrecision->30],",",N[h5[s],2],",",N[h4[s],2]]] 


k/j=+1 m=4 m=3 m=2 m=1 
-5 -9.49999986+0.20959825 i -9.49999986-0.20959825 i -9.90124704+2.1888866 i -10.90798882+2.2315190 i 
-4 -7.49999967+0.20959626 i -7.49999967-0.20959626 i -7.88383042+2.09028202 i -8.89331529+2.14210505 i 


3 
2 
“i 
es ee ee 


-5 
4 
3 
2 -4.49999794+0.29583878 i 
A 
0 ___| _-0.49930262-0.20834700i| _-0.49930262+0.20834701 | C-0.60999138-1.320191585 i] 


Table 11: Zeros of equation ¢(4, z) = 0 
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As we can see the program has 2 variables } = +/—1 and m = 1 + 4 and gives us 8 cases that have been 
listed in the table 11, in detail for 5 consecutive roots. Moreover, as we can see, there is no critical line for 
these roots (because Re(z) < 0), in the interval (0,1) but not in the interval greater than 1. Also the cases 
of transcendental equations ¢(2,z) = 0, ¢(3,z) = 0,¢(7,z) = 0 do not have the real part of complex zeros in 
the interval (0,1) have been checked. 


6.6.6. The zeros of equations ¢(8, z) = 0,¢(5,z) = 0,¢(3, z) = 0,¢(2, z) = O and the tables 
of these. 


Based on the same theory and method of chapter 6.5 we calculate the zeros around the interval (0, 1) 
below zero and above unity. If the values are not in this interval for the values of Re(z) then it means that 
there are no zeros in the interval (0,1) and therefore there are no critical lines other than 1/2. From the 
given examination we see that there are all points outside the interval and per pair of points of conjugate 
complexes. They have all in the interval (—1,0) from 1 pair upwards as q grows in the form ¢(q, z) that we 
consider. This means that as q increases in value, the probability of finding a critical line in the interval 
(0,1) is zero. The process of developing complex pairs is given on page 42 and they show a peculiarity 
around —1/2 and —1. 


I. Zeros of ¢(8, z) = 0. 


m=1 m=2 m=3 m=4 
es m=7 Pp mes 


-3.25350163-3.0535038027 i -3.25350163+3.053503802 i -3.49999972-0.74709613 i -3.49999972+0.74709613 i 
-0.87341271+2.4140275537 i -0.87341271-2.414027553 i -1.49998353-0.74710997 i -1.49998353+0.74710997 i 
7.040825660+8.2516499142 i} _7.040825660-8.2516499142 j 7.040825660-8.25164991 i 7.040825660+8.25164991 i 


-1, Se 0. <a i -1.50000043+0.3340980778 i 1. oe 0.0994660902i| 1.4999999681+0.0994660902 i 
1.3377*10*11-1.0538*10%11 i] _1.3377*10%11+1.0538*10%11i| _1.3777*10*11-3.1208*10*101| _1.3777*10%11+3.1208*10*10 1 
8.4870*10*11-1.2609*10%11i,} 8.4870*10%11+1.2609*10%11i,| _8.4693*10%11-3.9579*10%101| _8.4693*10*11+3.9579*10"10 | 


Table 12: Zeros of equation ¢(8, z) = 0 
II. Zeros of ¢(5, z) = 0. 
-2.49999557954879915011 | _-1.499979601926572792463 Tena 364843142 i en 364797934307 i nes 


-0.49977397564859132669 |  -0.49977397564859132669 | -0.500022339-0.3640594749i | -0.5003624239+0.36405780178i| __-0.685636087-1.6163407776 i 
707.606486559+1.625093353 i | 1.3780502*109+1.35785"1049 ||  424.34743499 -317.98956580 i 


pm me meme met 
“1,8318044283+1,87991492354| _ -0,685636087-1,61634077761 | -0.6856360874+1.61634077761| _-0.50036242-0.3640578017 -0.50036242+0,3640578017 i 

-0,685636087+1.61634077761| _22.85137203201030815180231| _*0.50000000+3,976353"10"8i | 1.4338*10%9-1.36540"10%9i | 1.43381602*10°9+1,36540083*109 
424,34743499 +317.98956580 | _424,34743499 -317.98956580; | _3,97635364883525332°10%8 i | _ 5.688577*10%9-1.6413*10"9i | _5.6885775*10%9+1,641333674*10"9 


Table 13: Zeros of equation ¢(5, z) = 0 


Al 


III. Zeros of ¢(3, z) = 0. 


Se Se << 
-2.4994443912584824 ~1,4987861940097745405 -2.63000774-1.267225286| -2.63000774+1.267225286 ~1,594811655534452208 ~1.49878619400823420555 


-0,4957156769130383 1710.309247510394492984 0,522828369-0,99303335 i -0,52282836+0,993033357 i 1710.3092475103944929 1710.3092475103944929 
1710.309247510394492984 1710.309247510394492984 1710.3092475103944929 1710.3092475103944929 1710.3092475103944929 1710.3092475103944929 


Table 14: Zeros of equation ¢(3, z) = 0 


IV. Zeros of ¢(2, z) = 0. 


~1.445569200+0.699260871065i | -1.4455692007+0.69926087106 i} _ -2.4599159118+0.7614749952i | _-2.4599159118-0.7614749952 i 


-1,445569200+0.69926087106 i | -1.4455692007+0.69926087106 i| -0.41213454795+0,5978119423 | -0.41213454795-0.5978119423 i 
71023*1095945.121569°10"S9 |_2.0740"10"60+9.7954*10'591 | 2.0740"10"60-2.7954%10"S91 | 


Table 15: Zeros of equation ¢(2, z) = 0 


Preceded by the given tables with zeros near the interval (0,1) laterally in order of magnitude. The 
colored zeros are those closest to the interval (0,1); notice immediately that they are divided into 3 cat- 
egories. Those where Re(z) is around —1/2 and those where Re(z) is close to —1. Following the logic of 
the matrices we construct the following Conclusion, which is general for the class of functional equations 
2(q,z), with gq >1 in Z+. 


6.6.7. Corollary 3. 


The zeros of the Generalized Riemann Zeta function ¢(q,z) =0,q >1 in Z+, do not be- 
long to the critical strip of the interval (0,1), but are outside. 


According to the Generalized root finding theorem [24] and the Infinite Periodic Radicals method, by 
which we solved the whole series of equations with g = 2,3,4,5,6,8 we observe that the roots, are divided 
into 3 categories, 2 in those where is around —1/2 and which approach it as q increases, and another one 
which moves away from —1/2 towards —1 as q increases. Of course we are looking at these roots which are 
their real parts in the interval (0, 1). 


We could also consider the case for gq in R in general, and such a brief examination will be done in section 
6.7 where we will consider the generalized equation G — ¢(z,q), i.e. we will also see some instances of ¢(q, z) 
with positive and negative gq. Thus we will have a complete picture of the functions ¢(q, z) and G — ¢(z,q) 
which are closely related to Riemann’s ¢(z). 


Let’s look in detail at the actual zeros per case in Table 16. 
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(3,2) 
-0.522828369+0.99303335 i] _-0.60999138+1.320191585i_| 
| C=C ~-0.4121345479-0.597811942 i, -0.522828369-0.99303335 i 
_—— = 4 
| Cd 
| 4 
[ie 


(8,2) 


-0.68563608+1.6163407776 i -0.7535529+1.8938854 i -0.87341271+2.4140275537 i 


-0.49977397564859132669 -0.500667-0.500214 i -0.500351168-0.747696558 i 
PY -0.4999612240.13404187 I -0.49999709+0.0994611864 i 
PY -0.49996122-0.13404187 I -0.49999709-0.0994611864 i 


-0.50036242+0.3640578017 i -0.7478539+1.8921442 i -0.500024879+0.334093564 i 


| sds -0.68563608-1.6163407776i | — -0.7535529-1.8938854 i -0.87341271-2.4140275537 i 
| «Ys -0,50036242-0.3640578017i | — -0.7478539-1.8921442 i -0.500024879-0.334093564 i 


-0.49977397564859132669 -0.500667+0.500214 i -0.500351168+0.747696558 i 


Table 16: Zeros table for g = {2,3, 4,5, 6, 8} 


Suppose we have the zeros 


1 12 2 11.22 22 fi 2) 3 1 .2 23 
neers } 5 {zie 12%, 422, } or) {aneizie 12h pees {zy 23s z3n } 5 
where {K = 1 +1, %1,%2,i3 € N} 
with ascending sequence of q, such that for the real parts of zeros apply the following 
Re | zi, )=—i-e! 
ak 27 “qk el, 0 Re( zt, } 3-2 
Gi, iy a 
2 
Re| 22, )=—1+€? Ee“, 0 
%, gt ak, 3S Sin Re ( 22, ) + -} 
eG > 0 tg 
Re[ 23, ) =-1+e63 3 
a, a, kK >> 2 Re ( z3, } > -1 
k,ne€Nw=l1lt+n ie 


Therefore all zeros they have a negative real part of the complex zero, which shows that no zero is in 
the critical Strip of the interval (0,1). So the Generalized Riemann Zeta equation ¢(q,z) = 0,q > 1 
in Z4, is outside of the critical line, but also has no critical line in the interval (0,1). For the Table 16, 
are used the functions Hurwitzzeta|q, z] or Zeta(q, z) to check the zeros with the program mathematica 12, 
clearly with more decimals to achieve approximations above 1071. 


The whole procedure for values of g > 1, clearly shows that the function ¢(q,z) = 0 does not generate 
zeros even just in the interval (0,1) which is the only critical line for the function ¢(z) = 0. Of course this 
does not exclude the possibility for other values of g, maybe it gives us some pairs or entire critical lines. 
This research will be completed below. 
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6.6.8. The zeros of equations ¢(—3,z) = 0,¢(—2,z) = 0, ¢(—2/3,z) = 0,¢(—1/2,z) = 0, 
¢(—1/4,z) = 0, and of equations ¢(1/4,z) = 0,¢(1/2,z) = 0,¢(2/3, z) = 0 and the tables of 
these. 


To be able to close the important issue of the zeros of ¢(q,z) = 0 and the positions of the values of 
Re(z) within the intervals (0,1), (—1,0)&(—oo, —1) we need to solve several forms of transcendental equa- 
tions and derive some characteristic matrices. The only method that has this possibility is the Periodic 
Radicals method, the same method that we solved the previous forms. the zeros we find will show the 
selective positions that are found and to which of these three intervals the Re(z) of each equation belongs. 


I. Control if Re(z) of zeros in interval (0, 1) 

I.1. Bernoulli polynomials 

Here there are 2 ways for integer values. Both the Infinity Periodic Radicals method and the Bernoulli 
polynomials. In this case we will prefer the Bernoulli polynomials. The Bernoulli polynomials are a gen- 


eralization of the Bernoulli numbers. They have a variety of interesting properties, and will feature in our 
proof of the Euler-Maclaurin Summation Formula. 


The generating function for the Bernoulli polynomials is the generating function for the Bernoulli numbers 
is multiplied by a term of e*¥. Our first observation of the Bernoulli polynomials is that the constant term 
of Bz(y) is in fact B. If we set y = 0, then B;,(0) = Bg. 


xery x 


em—1 er—1 


Proposition. The Bernoulli polynomials, Bz(y) satisfy the recurrence relation 


a>: (f ) Bout 


n=0 
3 By (y)a* xevy 
k! ez —1 
k=0 
x 
— ery 
et — 1 . 
Brak QS (ay)* 
= kl pe k! 
k=0 k=0 


We find the Cauchy product by the same process as earlier to obtain 


“8B ak oS zy)*-" Bra” 
a = ae 


Now, if we compare terms of x in the right and left-hand summations, we see the following: 


eS (5) Bout 


n=0 


as desired Using this recurrence, we can calculate the _rst few Bernoulli polynomials: 


Bo(y) =1 
Biy)=y-5 
r 1 
Bo(y) =y “iste 
3 
B = oe ed a 
3(y) =Yy a¥ + 5¥ 
1 
B = 49 3 Qi Se 
a(y) =Y Ue an 
5 5 1 
B = 5 4 3 
5(y) =Y ay + 5¥ — BY 


Notice our earlier result that the constant term of each polynomial is a Bernoulli number. Let us consider 
one of these polynomials, say Bs(y), more closely. What if we differentiate it, or integrate it over 0 to 1? 


d 1 
— Bs(y) = 5y* — 10y? + 5y? 
dy 6 
1 
= 5 4 i 9) 3 2 aaa te 
(u yo ry 30 
= 5Ba(y) 
so we observe that in this case, 
B,(y) = kBr-1(y) 
We also see , : 
1 5 1 
Bs(y) = sy® — sy? + sy" y? 
[ 6 ee 
=0 


These two observations are, remarkably, true in general. In fact, they give us the following inductive defi- 
nition for the Bernoulli numbers. 


More generally about Bernoulli polynomials apply: 


1. Boy) =1 
2. Bi(y) = kBr-1(y); 


3. fj Be(y)dy =0 for k > 1 
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I.2. Zeros of ¢(—3,z) = 0 and ¢(—2,z) = 0. 


In special functions, the following applies to the Zeta function 


=o Bq WN 1, z) 
¢( q, z) 7 q 5 1 
Therefore in these 2 cases we will have: 
B(3 
((-2,2) = - 2:3) 95 


= 1/3 (—(z/2) + (322) /2— 23) =0 


1 
the zeros are: z = O||z = giz =1 


And for the case ¢(—3, z) = 0 we have: 


(8,2) = 264) 


= 1/4 (1/30 — 2? + 228 — z*) =0 


the zeros are: 


1 2 1 2 
zg=-[1—4/1-2,/2]Jz2= = [144/1-2)/4] 

2 15 2 15 

1 | / 2 1 2 


z = 0.240335||z = 0.759665]| 
z = —0.157704||z = 1.1577 


=0> 


Solving other cases we come to the conclusion that we have 2 zeroes only real in the interval (0,1). Therefore 
we have no critical lines in the interval (0, 1). 


1.3. Zeros of ¢(—2/3, z) = 0, ¢(—1/2, z) = 0, ¢(—1/4,z) =0 


This case can only be solved with Infinity Periodic Radicals with use Generalized Theorem[24]. The matrix 
of the program that generates all solutions is the following: 
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The matrix program for ¢(—2/3, z) = 0,¢(—1/2, z) = 0,¢(—1/4, z) = 0 in mathematica with infinity 
Periodic Radicals is: 


Off[FindRoot: bddir.FindR oot::isine. Power::infy Infinity: :indet 


Genera 


h2[x_]:=Zeta[r.(x-1)]: 
h[x_]-=-1/(x-1))"*4Zetale.(«-1)] 
h3[x_]:=Zeta[r,x]: 
Forlk=-d,k<=d,k+, 
flu_]:=1+Exp[-((Log[u]+2 I zk)/r)]; 
xr=Nest{f[h2[#]]&,-1/100,1]; 
FQ2=N[xr, 10]: 
s=N[y1.40]/.FindRoot[h3[y1]= 0, {y1.FQ2}, WorkingPrecision->60]: 
If[Abs[h3[s]]<10’(-5), 
Print["x(",k,")=",s1=N[y1.40)/.FindRoot[h3[y1]= 0, {yl.s}. WorkingPrecision->60] 
-"."SN[h{s1]].".".N[h3[s1]]]:]] 


1. The zeros of ¢(—2/3, z) = 0 are: 


z, = 0.703226180818014 
Zq = 0.112695242215194 


2. The zeros of ¢(—1/2, z) = 0 is: 


z, = 0.066489124132138 


3. The zeros of ¢(—1/4,z) = 0 is: 


Z1 = 0.582245789118785 


As we observe in all three cases we have only real roots although they are located in the interval (0, 1) they 
obviously neither create nor belong to any critical line. 


1.4. Zeros of ¢(1/4, z) = 0, ¢(1/2, z) = 0, ¢(2/3, z) = 0 


Using the same program matrix we find all 3 zeros respectively. The matrix of the program that gen- 
erates all solutions is the following: 
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General::cvmit,ReplaceAll::reps. Divide 
d:=1:r:=-2/3 or -1/2 or -1/4; 
h2[x_]:=Zeta[r.(x-1)]: 
h[x_]=-1/((x-1))'r+Zetalr.(%-1)]; 
h3[x_]:=Zeta[r,x]: 

Forfk=-d.k<=d.k++, 
flu_]:=1+Exp[-((Log[u}+2 I zk)/r)]: 
xr=Nest[f[h2[#]]&.-1/100,1]: 
FQ2=N[xr, 10]: 


s=N[y1.40]/.FindRoot[h3[y1]= 0, {y1.FQ2}, WorkingPrecision->60]: 


If[Abs[h3[s]]<10°(-5), 


Print["x(".k.")}=".s1=N[y1.40)’.FindRoot[h3[y1]= 0. {yl.s}. WorkingPrecision->60] 


."S"SN[h[s1]].".".N[h3[s1]]]:]] 


1. The zeros of ¢(1/4,z) = 0 are: 


z1 = 0.408064210108089 


2. The zeros of ¢(1/2, z) = 0 is: 


21 = 0.802721828598366 


3. The zeros of ¢(2/3,z) = 0 is: 


Z, = 0.221547787091544 


Therefore for all cases with q < 0 and q in (0,1) while there are zeros in the set of real numbers in the 


interval (0,1) there are no critical lines in this interval. 


6.6.9. Dirichlet L-Series and their relationship to the critical line. 


I. Dirichlet L-series is a series of the form [27] 


Lx(s,x) = >> Xe(n) 07° 
n=0 


where the number theoretic character x(n) is an integer function with period k, are called Dirichlet L-series. 
These series are very important in additive number theory (they were used, for instance, to prove Dirichlet’s 
theorem), and have a close connection with modular forms. Dirichlet L-series can be written as sums of 
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Lerch transcendents with z a power of e?**/*, Dirichlet L-series is implemented in the Wolfram Language as 
DirichletL [k, 7,5] for the Dirichlet character y(n) with modulus & and index j. The generalized Riemann 
hypothesis conjectures that neither the Riemann zeta function nor any Dirichlet L-series has a zero with 
real part larger than 1/2. 


The Dirichlet lambda function also is, 


ie ui) Care) 


The function Dirichlet (k,j,s), gives the Dirichlet L-function L(y,s) for the Dirichlet character y(n) with 
modulus /& and index 7. 


Example(with mathematica): 


FindRoot [Dirichlet [8, 3, t] == 0, {t, 14}], Find Root|[Dirichlet([8, 3, t] == 0, {t, 15/}] 
Out : {t + 0.5 + 12.9881i}, {t + 0.5 + 16.3426i} 


We therefore always have as a critical line 1/2. 


II. The Dirichlet beta function is defined by the sum 


= ~ pals _ 9-8 1 
60) = 2 Gagiap ~? »(-1,55) 


where ®(z,s,a) is the Lerch transcendent. The beta function can be written in terms of the Hurwitz zeta 


function ¢(x, a) by 
woe F(a) <0) 


The beta function can be defined over the whole complex plane using analytic continuation, 


B(1—z) = (2) sin (5-7-2) -P(z) - B(z) 


where I'(z) is the gamma function. 


Example(with mathematica): 


Find Root[DirichletLBeta[s] == 0, {t, 147}] 
Out : {t > 0.5 + 12.98817} 


We therefore always have as a critical line 1/2. 
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III. Dirichlet Eta Function 


The Dirichlet eta function is the function 7(s) defined [28] by 


nfs) = = (1-289) 669) 


n=1 


where ¢(s) is the Riemann zeta function. Note that Borwein and Borwein (1987, p. 289) use the nota- 
tion a(s) instead of n(s). The function is also known as the alternating zeta function and denoted ¢*(s), 
n(0) = 1/2 is defined by setting s = 0 in the right-hand side of (2), while 7(1) = In2 (sometimes called 
the alternating harmonic series) is defined using the left-hand side. The function vanishes at each zero of 
1 — 2!~* except s = 1. The eta function is related to the Riemann zeta function and Dirichlet lambda 
function by 


2 ot 2 
And ¢(v) + 7(v) = 2A(v). The eta function is also a special case of the polylogarithm function, 


n(a) = — Lie (—1) 


The derivative of the eta function is given by 


nf (x) = 2'-*(In2)¢(x) + (1 — 2°*) C'(a) 


Therefore, the natural logarithm of 2 is 


Example(with mathematica): 


FindRoot [DirichletLEeta [s] = =0,{t,14I}], 
Out : {t > 0.5 + 14.1347i} 


We therefore always have as a critical line 1/2. 


IV. Ramanujan z-Diriehlet series 


In this paper we consider the zeros of the Dirichlet series 


(c = Res > 13/2) 


» T(n)z” =z PsP (2215 


It is known that L(s) has all of its non-real zeros in the strip 11/2 < 0 < 13/2 and it is conjectured that all 
these zeros are on o = 6 and are simple. It has been shown (by Hafner [29]) that a positive proportion of 
the zeros of L(s) are of odd multiplicity and are on o = 6, but it is unknown whether any of the zeros are 
simple. Some other Dirichlet series are known to have infinitely many simple zeros. Heath-Brown [80] and 
Selberg independently observed that the work [31] of Levinson in 1974 implies that a positive proportion of 
the zeros of the Riemann zeta function are simple. Of course, the same result works for a Dirichlet L-function. 


Example(with mathematica): 


FindRoot[RamanujanTauL|s]==0,{t,5+9I}], 
Out : {t > 6 + 9.222381} 


We therefore always have as a critical line in 6. 


V. Zeros of the Davenport-Heilbronn counterexample 
We compute zeros off the critical line of a Dirichlet series considered by H. Davenport and H. Heil- 
bronn. This computation is accomplished by deforming a Dirichlet series with a set of known zeros into the 


Davenport- Heilbronn series. 


For s=a+it with o > 1, let 


bates (1) 


be a Dirichlet series with periodic coefficients of period 5. Then f;(s) defines an entire function satisfying 
the following functional equation 


f(s) =T~**? x0(s)f(1— 8) (2) 


with T = 5 and 


x2(s) = 2(2m)82P(1 — 8) cos (=) (3) 


In 1936, H. Davenport and H. Heilbronn (see [32]) showed that f1(s), as defined in (1), has zeros off the 
critical line 0 = 1/2. In 1994, R. Spira (see [33]) computed the following zeros of the Davenport-Heilbronn 
example: 


0.808517 + 85.6993482, 0.650830 + 114.1633432 
0.574356 + 166.4793062, 0.724258 + 176.7024612 


In this note we present a scheme for computing additional zeros of the Davenport- Heilbronn Dirichlet 
series. We certainly see isolated cases of pairs of complexes that individually define critical lines. We have 
seen this before in other cases with ¢-functions and corresponding functional equations. 


”It is not a rare phenomenon but here we have to clarify which equations the Riemann 
hypothesis accepts as real cases as true cases that define the structure of the given hypothe- 


sis” 
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6.6.10. Solve the generalized equation G — ¢(z,q) = 0 of Riemann. 


To get a general picture of all forms of ¢() functions we have to deal with the general case which we 
will call G — ¢(z,q). If we define a complex z = o + Ti and D is the set of real o, it will converge absolutely 
for 0 > 1. For the value z = 1 it will diverge. So it will have 2 poles z = 1 and where the denominator is 
zeroed. The form of the sum will now be, 


oS 1 a+ 8 
= = wh = wl 
G— ¢(z,q) », mee ,{a,8 €R},2€C (G.1) 
Therefore we have 2 poles, z = 1 and z = —8/a. There are 4 forms ( of functions that are equivalent 


to the generalized function of G — ¢(z,q), which means, that they are not always equal. The last 2 provide 
the solution of the generalized equation G — ¢ (z,q) = 0. The Important point of solving G — ¢(z,q) = 0 is 
to check if there are zeros with 0 < 0 < 1, i.e. if Re(z) is within the open interval (0, 1). 


Equivalent relations functions that hold with respect to the values a,b and q. 


If where g = ate ,{a, 8 € R},z €C, apply per case, 
G—¢(z,q) =a *- HurwitzZeta(z,q) (1) 
G—C(z,q) =a7*- Zeta(z,q) (2) (G.2) 


G-—C(z,q)=a°*- ((-7y) + Zeta(z,q— 0) = T, — Zeta(z,q) (8) 


G—¢(z,q)=a*- (2) + Zeta(z,q+ ) = Ty — Zeta(z,q) (4) 


But for the process of solving with the Infinite Periodic Radicals method and the generalized G.Rt theorem 
we need to analyze further the cases related to the values of a, @ and for which functions are equal to the 
G—(¢(z,q). The main thing we are mainly interested in are the functions T; — Zeta (z, q) & T2 — Zeta(z, q) 
because they are represented by exponential functions with respect to q and therefore give multiple different 
zeros. In total we distinguish 5 cases related to the signs of a, @ and q. 
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I. The case where a < 0,3 >Oandq>0 & |q| <1. 


In this case the following equations apply: 
1 z 
G— ¢(z,q) =a * - HurwitzZeta(z,q) =a-*- ((<) + Zeta(z,q+ 1) = To — Zeta(z,q) VzEC (G.2) 


Since the objective is G— ¢(z,q) = 0 with Ty — Zeta(z, q) = 0 implies: 


G—C¢(z,q) =a~*- HurwitzZeta (z,q) =a*- (2) + Zeta(z,q+ ) =0 
with Generator (+ ae => 

enerato (4) u face 
{k € Z,u 40,q 4 0} and 
= _2-i-m-k | log(u) 
2= I) ~ Toate) * Tos/a) 


, Generator function. 


Program(example) G.¢.1. The matrix program in mathematica with infinity Periodic Radicals is: 


Off[FindRoot::bddir, FindRoot::jsing, Power::infy, Infinity::indet, 
General::cvmit, ReplaceAll::reps, Divide::infy, FindRoot::Istol,N::meprec]; 
d:= 10; a :=-3; b:=2; q :=(at b)/a; 

h2[x_] := -Zeta[x, (1 + q)]; 

h[x_] = (1((q))x + Zeta[x, (1 + q)]; 

h3[x_] := HurwitzZeta[x, q];j = 0; 

For[k = -d, k <= d, k++, 

f[u_] = (2 I\[Pi] k)/Log[1/q] + Log[u]/Log[1/q]; 

xr = Nest[f[h2[#]] &, -1/1000, 2]; 

FQ2 = N[xr, 10]; 


s=N[yl, 40] /. 

FindRoot[h[y1] == 0, {yl, FQ2}, WorkingPrecision -> 180]; 

If]0 < Re[s] <1, 

j=jt1; Print["x(", k, "=", sl =N[y1, 40] /. FindRoot[h[y1] = 0, {y1, s}, WorkingPrecision -> 180], ",", 
N[h[s1]],",", 


N{h3[Conjugate[s1]]], ",", N[h3[s1]]];]] ; 
Print["Percent of Zeros %=", N[j/(2 d+ 1)*100]] 


This case is very basic like the next one because it produces zeros belonging to the interval (0,1) that are 
doubly conjugate and in particular 0 < Re(z) < 1/2. They indicate that the hypothesis is rejected if its 
formulation refers to the most generalized form of the z-function. Of course, this fact if it holds depends on 
the assumptions we make in the definition of the Hypothesis from the beginning. 


Some results of the project are given in the table 17: 
Example a = —3 & 8 = 2. 
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number Zeros 


0.38273620158054312 + 59.88327733641703154 i 
0.38273620158054312 - 59.88327733641703154 i 
0.14082903205802571 + 53.3070632589853913 i 
0.14082903205802571 - 53.3070632589853913 i 
0.09375600099886612 + 46.080744571745147 8 i 
0.09375600099886612 - 46.080744571745147 8 i 
0.47265057875792859 + 39.6960570839325771 i 
0.47265057875792859 - 39.6960570839325771 i 
0.46209119678156450 + 33.64347641517712960 i 
0.46209119678156450 - 33.64347641517712960 i 


G-U(z,q)=0, {a<0,B>0,q>0, |q|<1} 


Table 17: 5 zeros for Re(z) in interval (0, 1) 


Therefore, as we observe, we have a number of conjugate zeros in the interval (0,1). This means that 
we will have scattered critical lines in pairs of complex zeros, which are infinite. The percentages cording 
to the number of zeros requested for more than d = 10 ranges from 76% to 99.99...% as the number of zeros 


requested increases. 


II. The case where a < 0,38 <0 andq>0 & |q| >1. 


In this case equations apply: 


G—¢(z,q) =a*- HurwitzZeta(z,q) = a7 


G—¢(z,q) =a*- Zeta(z,q) Vz EC 


G—C(z,q) =a? -HurwitzZeta(z,q) = a7?-((4 ” + Zeta z,q+1)) = Te —- Zeta(z,q 
q 


Zz, (- (4) + Zeta(z,q+ ) = T — Zeta(z, q) 


Since the objective is G — ¢(z,g) = 0 with T2 — 


Zeta (z2q) = 0 implies: 


z 
with Generator (3) =u> 


{k € Z,u 40,q 4 0} and 
fos log(u) 


G—C(z,q) =a-*-+ HurwitzZeta (z,q) =a-~*-((4 "£768 z,q+1)) =0 
q 


log(1/q) _ log(1/q) 


, Generator function. 


Program(example) G.¢.2. The matrix program in mathematica with infinity Periodic Radicals is: 


Some results of the project are given in the table 18: 


Example a=—-3 & B=-1 
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Off[FindRoot::bddir,FindRoot::jsing,Power::infy Infinity::indet,General::cvmit,ReplaceAll::reps,Divide::inf 
y,FindRoot::Istol,N::meprec]; 
d:=20;a:=-3;b:=-1;q:=(a+b)/a; 
h2[x_]:=-Zeta[x,(1+q)]; 
h{x_]:=(1/(q))x+Zeta[x,(1+q)]: 
h3[x_]:=HurwitzZeta[x,q]:j:=0; 
For[k=-d,k<=d,k++, 
flu_]:=(2 Im k)/Log[1/q]+Log[u]/Log{[1/q]; 
xr=Nest[f[h2[#]]&,-1/1000,2]: 
FQ2=N[xr, 10]; 
s=N[y1,40]/.FindRoot[h[y1]==0, {y1,.FQ2},WorkingPrecision->180]; 
If[0<Re[s]<1 j=j+1 ;Print["x(",k,")=",s1=N[y1,40]/.FindRoot[h[y1]==0, {y1,s}, WorkingPrecision- 
>180],",",N[h[s1]],",",N[h3[Conjugate[s1]]],",",N[h3[s1]]];]]; 
Print|"Percent of Zeros %=",N[j/(2 d+1)*100]] 


| SC _0.8200108859434798 - 17.2889024511240770206 i_| 
| (0.373 7554272435496 - 24.8603622893908232381 i_| 


| | 0.3902145457078484 - 46,3444286577441579380 i_| 
| 4 —(0.4004936636741982 + 62.8931327281430963430i_| 
| C0. 4004936636741982 - 62.8931327281430963430i_| 


0.61684575102099114 - 78.307731086389698349 i 


Table 18: 5 zeros for Re(z) in interval (0, 1) 


So here too we have a number of conjugate zeros in the interval (0, 1).This means that we will have scattered 
critical lines in pairs of complex zeros, which are infinite. The percentages from d = 20 range from 40% to 
99.99...% as the number of required zeros increases. 


III. The case where a > 0, 8 < 0, |a| > |G| andqg >0< & |q| < 1 & The case wherea > 0,8 < 0, 
|a| < |B] and gq >0< |q| <1 


They have no solution because the G — ¢ (z, g) tends to infinity: 


= (a-k+ 8) 


= 1 
G-¢(z,q)=5 7 + 00, where q= 8 {a,8 € R},z C 
k=1 


We therefore do not examine these cases. 
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IV. The case where a > 0, 8 >O and q>0& |q|>1 


In this case equations apply: 


G-—¢(z,q) =a *- HurwitzZeta(z,q) =a~*- (2) + Zeta(z,q+ ) = T2 — Zeta(z, q) 
G-—¢(z,q) =a *- HurwitzZeta(z,q) =a~*- (- (4) + Zeta(z,q+ ») = T, — Zeta(z, q) 


q-l 


G—¢(z,q) =a*-Zeta(z,q) VWzeEC 


Since the objective is G — ¢(z,q) = 0 with T2 — Zeta(z,q) = 0 implies: 


G —C(z,q) =a-*- HurwitzZeta (z,q) =a~*- (2) + Zeta(z,q+ \) =0 


z 
with Generator (4) =u> 


{k € Z,u#0,q 4 0} and 


z=f(u= esti yay + AB Generator function. 


Program(example) G.¢.4. The matrix program in mathematica with infinity Periodic Radicals is: 


Off[FindRoot::bdd ( ising. Power: :infy. Infinity: :indet 
A bet haere 
General::cvmit.ReplaceAll::reps.Divide: :infy.FindRoot::lstol.N::meprec.Sum::div]: 


d:=35:a:=3:b:=7:p:=(atb)/a: 
h2[x_]:=-Zeta[x.(1+p)]: 
h[x_]:=(1/(p))*x+Zeta[x.(1+p)]; 
h3[x_]:=HurwitzZeta[x.p]: 


= 1 : 
= 2) > J; 
nae = 2 Gap 


Forfk=-d.k<=d.k++, 
fju_J=( 1 7k)/Log[1/p]}+Log[u}/Log[1/p]: 
xi=Nest[f[h2 [#]]&.-10%(-40).1]: 
FQ2=N[xr.10]: 
s=N[y1.40]/.FindRoot[h4 [yl ]== 0. {y1.FQ2}.WorkingPrecision->5 80]: 
If{Abs[h4[s]]<10“(-40) j=j+1: 
Print["x(".k,")}=",s1=N[y1,40]/.FindRoot[h4 [yl ]== 0, {y1.s}. 
WorkingPrecision->680] .".", N[h[s1]].".".N[h4[s1]].".".N[h3[s1]]]:]] : 
Print["Pososto %=".N[j/(2d+1)*100]] 


Some results of the project are given in the table 19 
Examplea=3 & B=7 
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(G.7) 


(G.8) 


G-U(z,q)=0, {a>0,B>0,q>0, |q|>1} 
1 ee 


85.47509383 


86.18462479+10.42600498 | 
86.18462479-10.42600498 | 


Ee ee 85.97351095+15.35412591 | 


Table 19: 5 zeros for Re(z) in interval (0, 1) 


Note that the zeros all with Re(z) outside the interval (0,1) per pair. We also notice that the values 
of Re(z) are very close and there is also a net true root at 85.475..about. 


V. The case where a < 0, 8>Oandq <0 & |q|>1 


The equations that apply in this case are: 


1 Z 
G—¢(z,q) =a”? - HurwitzZeta(z,q) = a~* - ((<) + Zeta(z, q+ 1) = T2 — Zeta(z, q) 


1 * : 
G — ¢(z,q) =a * - HurwitzZeta(z,q) =a™*- (- (=) + Zeta(z,q+ 0) = T, — Zeta(z, q) (G2) 
q= 
G— ¢(z,q) =a? - Zeta(z,q) VWzEC 
Since the objective is G — ¢(z,q) = 0 with Tz — Zeta(z, q) = 0 implies: 
= ‘ _ Ive 
G - ¢(z,q) =a * - HurwitzZeta(z,q) =a-* - ((<) + Zeta(z,q+ 1) =0 
ith Generato Sy, > 
wil ner 1 ea (oes =U 
q (G.10) 


{k € Z,u40,q 4 0} and 
a (es log(u) 
ee INO Tos(/q) | log(1/a) 


, Generator function. 


This is a peculiar case we are considering, which has no double conjugates but only simple conjugates, 
so simple points. It is not a rejection of Riemann’s Hypothesis because no critical lines are defined. Obvi- 
ously we can look at many more examples and draw some conclusion about the region we are considering. 
But for cases I&II we have drawn very useful conclusions that warn us that the Riemann Hypothesis bends 
when it faces the generalized function G — ¢(z,q). 
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Program(example) G.¢.5. The matrix program in mathematica with in_nity Periodic Radicals is: 


Off[FindRoot::bddir FindR oot::jsing, Power::infy Infinity::indet,General::cvmit ReplaceAll::reps Divide::inf 
y.FindRoot::Istol,N::meprec]; 

d:=20;a:=-3:b:=10:p:=(atb)/a: 

h2[x_]:=-HurwitzZeta[x.(1+p)}; 

h{x_]=(1/(p))’x+HurwitzZeta[x.(1+p)]: 

h3[x_]:=a’(-x)*HurwitzZeta[x,p]:j:=0: 


x 


1 
bis I= ape 


k=! 
Forlk=-d.k<<d,k++.flu_]:=(211k+Log[u])/Log[1/p]: 
xr=Nest[f[h2[#]]&.1/1000,1]: 
FQ2=N{xr, 10]; 
s=N[y1,40]/.FindRoot[h4[y1 F=0,{y1,FQ2}.WorkingPrecision->180]: 
If[0<Re[s]<1,j=j+1; Print["x(",k,")=",sl=N[y2,40]/ FindRoot[h4[y2 =0,{y2,s}, 
WorkingPrecision->180],".".N[h{s1]].".".N[h4[s1]}.".".N[h3[s1]]]:]:] 
Print|"Percent of Zeros %=",N[j/(2 d+1)*100]] 


Some results of the project are given in the table 20 
Example a = —3 & a= 10 


G-4(2,q}=0, {a<0,B>0,9<0, |q|>1 
0.4366208398066607628 + 0.47787514549539380431 i 
| =.=: 0,0151755073772275510 - 7.53810974894711608176 i 


| 3. 0,419986170878482310452 - 11.315286315927897635 | 
| = 4_~—_| 0,3279682235408925957 - 1.918505421228014939493 i 


Table 20: 5 zeros for Re(z) in interval (0, 1) 


We see that the zeros of this class are all with Re(z) within the interval (0,1) per root. We also note that 
the values of Im(z) do not have opposite signs are what we call orphan complex roots it is now obvious that 
for the generalized equation G — ¢(z,q) = 0 its zeros where R(z) is in the interval (0,1), correspond only to 
cases I,]I,V. But this equation is derived from the generalized z() function, so the Riemann hypothesis is 
Rejected in these cases as well. 
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Epilogue to Chapter 6 


I. The zeros of ¢(z) = 0 


The fact that the roots of ¢(z) = 0 all lie on the critical line, follows from 5 proofs, both basic and el- 
ementary. 


1.Using direct method of inverse of ¢() on functional equations, chapter 6.1.1. 
2.Solutions of the Equation ¢(s)/¢(1 — s) = 1, s is Complex number, chapter 6.1.3. 
3.With general inverse random function theory, chapter 6.2. 

4.Statistical evidence method chapter 6.3. 

5.With the define the function A(N), chapter 6.4.1. 


Therefore, the proof that all zeros of ¢(z) = 0 belong on the critical line is complete, since the func- 
tion ¢(z) = 0 as proved many times has zeros only on the critical line which means that Re(z) = 1/2 and 
only this one. 


II. Other variations of the function ¢(z) = 0 


Finally, the cases ¢(z,q) = 0 or ¢(q,z) = 0 or G — ¢(z,q) = 0, q € R, are considered in detail in these 4 
intervals to clarify whether and where there are zeros and whether they correspond to critical lines. The 
analytical examination yields zeros with Re(z) in the interval (0,1) for the cases ¢ (q * z) = 0 where gq > 1/2 
in R, G— ¢(z,q) = 0 for at least 2 general cases out of 5, and for case Davenport Heilbronn case, which 
has zeros of a Dirichlet series periodic. Therefore, whether or not the Riemann Hypothesis holds depends 
on whether or not we accept the fact if we consider only ¢(z) = 0. If in the definition of the Hypothesis 
the other forms or part of them also hold i.e. ¢ (q¢ * z) = 0 or ¢(q,z) = 0 or G- C(z,q) = 0 with q E Rin 
general or the Epstein Function or the Dirichlet series and others, then the Riemann Hypothesis does not 
hold in general, but only in special cases which are stated in detail. But we have to accept an indisputable 
fact that in order to have a visual perception of which ¢-functions have zeros with Re(z) = 1/2 or have 
Re(z) different from 1/2 or other critical lines outside the interval (0,1), we have to solve transcendental 
equations [24] with accuracy. 


III. In which cases the Riemann case is valid or rejected. 


Cases I, II of the characteristic generalized equation G — ¢(z,q) = 0,¢ > 0,|¢| < lor |¢] > lade R 
(pages 52-58) as well as the case ¢ (q* z) = 0 with g > 1/2 (pages 30-32) are subclass of the generalized 
one mentioned above, and some cases of the Dirichlet periodic series Davenport-Heilbronn (page 51) reject 
the Riemann Hypothesis. So if the hypothesis is formulated as follows: ”For any ¢-function, the zeros of 
the resulting equation coincide on the critical line with Re(z) = 1/2, the formulation is complicated as a 
Latent Hypothesis and does NOT hold”. But if it is formulated as ” The zeros only of ¢(z) = 0 coincide on 
the critical line for Re(z) = 1/2 is correct and the Hypothesis is Valid”, as has been shown. Obviously this 
happened because the equations of other related functions i.e other related equations arising from related 
functions with the ¢-function were not completely or adequately or at all and solved and there was a large 
knowledge gap related to where their zeros are located with respect mainly to the interval (0,1). These 
equations are forms of Transcendental equations quite difficult and to solve them a strong method and a 
generalized Theorem for solving such equations is needed. 
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Part II. Distribution of primes numbers and Zeros 
7. Introduction 


One way to get a feeling for where the primes are within the natural numbers is to count the number 
of primes a(n) less than or equal to a given number n. There is no real hope for finding a simple rule 
for t(n), but the Prime Number Theorem1 (PNT) says that a(n) can be approximated by the logarithmic 
integral function. The distribution of prime numbers is most simply expressed as the (discontinuous) step 
function 7(x), where 7(a) (Fig.9), is the number of primes less than or equal to x. 


Figure 9: The primes of function 7() 


The function 7(a) in relation to the random int. x. It turns out that (a) can be expressed exactly as the 
limit of a sequence of smooth functions R,,(a). To define R,,(x) we first introduce the logarithmic integral 
function Li(2), 


which appears throughout the analytical theory of the prime distribution: 


This is a smooth function which simply gives the area under the curve of the function 1/logu in the 
interval [2,2]. Don Zagier explains the reasoning behind the function Li in his excellent introductory ar- 
ticle “The first 50 million prime numbers” [14], based on his inaugural lecture held at Bonn University, 
May 5, 1975: “A good approximation to a(a)”, which was first given by Gauss is obtained by taking as 
starting point the empirical fact that the frequency of prime numbers near a very large number z is almost 
exactly 1/loga. From this, the number of prime numbers up to x should be approximately given by the 
logarithmic sum Li(x) = 1/log2 + 1/log3 + 1/logz or, what is essentially the same, by the logarithmic 


integral. 
“ d 
Li(e) = f 
2g logu 


Using Li(a) we then define another smooth function, R(x), first introduced by Riemann in his original eight 
page paper, and given by 


Riemann’s research on prime numbers suggests that the probability for a large number « to be prime should 


be even closer to 1/log a if one counted not only the prime numbers but also the powers of primes, count- 
ing the square of a prime as half a prime, the cube of a prime as a third, etc. This leads to the approximation: 


if 1/2 1 1 ~~ : 
m(a) + 5n(a P) + a(a 3) 4. = Li(a) 


or, equivalently [by means of the Mobius inversion formula] 


(a) © Lila) — 5 Li(a'/?) — S Lia’) - 


The function on the right side of this formula is denoted by R(x), in honour of Riemann. It represents 
an amazingly good approximation to 7(x). For those in the audience who know a little function theory, 
perhaps I might add that R(x) is an entire function of log x, given by the rapidly converging power series: 


1430 (In x)* 
« REIC(E + 1) 
where ¢(& + 1) is the Riemann zeta. 
pen (Ronkaivin. Zéa9s) 
taey v2 + (87.58. 
$54 12 92.93..9/ 
307 $F V2 + (30.42.97 
$4 12+ 5.01 yf 


20% t= 4. 4/2 + (21.02...)f 
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ZF 
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Figure 10: The non-trivial zeros in the Critical Line 


Here we see the zeros of the Riemann zeta function in the complex plane. These fall into two categories, triv- 
ial and nontrivialzeros (Fig.10). Here are some tables on nontrivial zeros compliled by Andrew Odlyzko[4]. 
The trivial zeros are simply the negative even integers. The nontrivial zeros are known to all lie in the 
critical strip that is 0 < Re[s] < 1, and always come in complex conjugate pairs. All known nontrivial 
zeros lie on the critical line Re[s] = 1/2. The Riemann Hypothesis states that they all lie on this line. 
The difference between the prime counting function and its “amazingly good approximation” R(x), i.e. the 
fluctuations in the distribution of primes, can be expressed in terms of the entire set of zeros of zeta, which 
we shall represent by p, via the function R itself: 


R(x) —1(x) = 5° R(x*) 


Obviously some of the x? are complex values, so here R is the analytic continuation of the real-valued 
function R defined previously. This was mentioned above by Zagier, and is known as the Gram Series 
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expansion: 
co 


VRC HR) 


k=1 
The contributions from the complex-conjugate pairs py and p_~, = px cancel each others’ imaginary parts, 
so 


n(x) = R(x) — S> R(x?) + N° Ti (a) 
k=1 


m=1 


Where T;,(x) = —R(ax?*) — R(ax?-*) are real-valued. We can now define the sequence of functions R,,(x) 
which approximate 7(2) in limit: 


n 


R,(x) = R(x) - S- R(a~?™) + S 0 T(x) (1) 


m=1 k=1 


f(x) = —2Re (>. Et (pr ro) + es i 5 loz(b) dt — log(2) (2) 


k=1 


where p, is the x*® complex zero of the zeta function. In this formula, Ei(z) (Mathematica’s built-in 


function ExpIntegral Ei(z) is the generalization of the logarithmic integral to complex numbers. These 
equations come from references [13], [14], and [15]. First, let M be the smallest integer such that 2!// < 2. 
We need to add only the first M-1 terms (that is, n = 1,2,..,/) in the sum in equation (1). For each of 
these values of N, we use equation (2) to compute the value of f (al/ ie However, we will add only the first 
N terms (that is, « = 1,2,..,. N) in the sum in equation (2). Because the purpose of this Demonstration is to 
show how the jumps in the step function m(a) can be closely approximated by adding to R(«) a correction 
term that involves zeta zeros, we ignore the integral and the log2 in second equation; this speeds up the 
computation and will not noticeably affect the graphs, especially for « more than about 5. The more zeros 
we use, the closer we can approximate m(). For larger x, the correction term must include more zeros in 
order to accurately approximate 7(). 


8. The GUE hypothesis. 


While many attempts to prove the RH had been made, a few amount of work has been devoted to the 
study of the distribution of zeros of the Zeta function. A major step has been done toward a detailed 
study of the distribution of zeros of the Zeta function by Hugh Montgomery [18], with the Montgomery pair 
correlation conjecture. Expressed in terms of the normalized spacing. 


log (Yn/(27)) 
2 


On — (Yn41 Yn) 


this conjecture is that, for M — oo 


b b . 2 
qin N+1snsN4+M5y € (as) [ p(O,u)du~ f 1- (=) a 


TU 


Tu 
It was first noted by the Freeman Dyson, a quantum physicist, during a now-legendary short teatime ex- 


change with Hugh Montgomery[18], that this is precisely the pair correlation function of eigenvalues of 
random hermitian matrices with independent normal distribution in (Fig.11, 12) of its coefficients. Such 
random hermitian matrices are called the Gauss unitary ensemble (GUE). As referred by Odlyzko in [16] for 
example, this motivates the GUE hypothesis which is the conjecture that the distribution of the normalized 
spacing between zeros of the Zeta function is asymptotically equal to the distribution of the GUE eigen val- 
ues. Where p(0, wu) is a certain probability density function, quite complicated to obtain for an expression of 


; i : sin 7u 
In other words, the density of normalized spacing between non-necessarily consecutive zeros is 1— ) i 
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2 4 
it). As reported by Odlyzko in [16], we have the Taylor expansion around zero. p(0,u) = ow - 2Tul +... 
which under the GUE hypothesis entails that the proportion of 6 n less than a given small value 6 is asymp- 


2 
TT 
totic to Go + O (5°). Thus very close pair of zeros are rare. 


Density 


Normalized spacing 


Figure 11: Probability density of the normalized spacing 


between non-necessarily consecutive zeros and the GUE prediction[20] 


— 


Figure 12: Probability density of the normalized spacing dn and the GUE prediction[20]. 


9. Gaps between zeros 


The table below lists the minimum and maximal values of normalized spacing between zeros 6n and of 

On + On41, and compares this with what is expected under the GUE hypothesis. It can be proved that 
2 4 

p(0,t) have the following Taylor expansion around 0, p(0,u) = 1 we & 


3 U 25 u* +--+ so in particular, for 
small delta 


6 2 
Prob(6n <6) = p(0,u)du ~ =o 
0 


so that the probability that the smallest dn are less than 6 for M consecutive values of dn is about. 


nr? a 1? 
1= (1 = = ~1-—exp (-F2%m) 


This was the value used in the sixth column of the table. The result can be also obtained for the 6,+6n41.-.. 


6 
Prob(5n + Sn41 <6) ~ se 
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from which we deduce the value of the last column (Table 21). 


Height Mini 4,, n | Mini Maxi Prob min | Prob min 6, + 
y j J 5 6, in GUE | 6,4; in GUE 


0.0005330 
0.0009764 
0.0005171 
0.0005202 
0.0006583 0.0966 
0.0004390 0.1080 
0.0004969 t 0.0874 
0.0004351 4.268 0.1067 
0.0004934 316 0.1019 
0.0008161 . 0.1060 
0.0004249 4.30¢ 0.1112 
0.0002799 15 0.0877 


0.1097 
0.1213 
0.1003 
0.1029 


Sash 


fs be ho et OF 
COSC ) 
-~OnN 


Table 21: The values of bn + 6n41 


For very large spacing in the GUE, as reported by Odlyzko in [16], des Cloizeaux and Mehta [17] have 
proved that 
log p(0,t) ~ —17t?/8 (t + co) 


which suggests that 
1/2 
ie i ee 
N+1<n<N+M 1 


Statistics of False zetazeros, 6-intervals and count of primes. 


10. General equations. 


With this statistic we find the crowd of individual parts intervals defined by consecutive Zetazeros in a 
fixed integer interval. Here we use 6 = 1000. We have the general Equation 6+ Zetazero (kin) - Ze- 
tazero (kf) = 0 and after given initial value in kin, i calculate the kf usually by the Newton method. 
From the initial and final value of kin, kf and by performing the process of successive intervals, i calcu- 
late the number(count) of the primes ones that are within the intervals. In this way will have consecutive 
intervals. S;, = { Zetazero (k + 1) < p < Zetazero (k),k € (kin, kf)} from where resulting the number of 
the primes in the given consecutive interval. The sum of the primes in the given interval 6 = 1000 will be 
obvious. 


ky 
an eer 


i=kin 


The False intervals will be in a normalized form Fs = 1000 — Ss over N >= 4000. 
11. The statistics 


The First Statistic looking for ways to show us which function is the most ideal to get closer the points of 
interest, uses the NonlinearModelFit method of the function y = (a+ b- )/(log(d + c- a) and therefore 
after determining the variables {a, b, c,d} we are able to make statistical and probable prediction at higher 
levels of numbers . This function is directly related to the function 7(x) = a/log« which was reported in 
the introduction that is, defining the number of primes numbers relative to «. We will do a double statistic 
of the intervals {6, False intervals } and {6, number of primes } that we are ultimately interested in 
the statistics mainly the count of the primes inside at 6 -intervals. 
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11.1 The first statistic is about count of False intervals and the count of the primes that correspond- 
ing in them. The table below gives it in aggregate until the count of primes equal ps = 29 corresponding to 
the count of False Fs = 971 and the Total range Integers for a interval D = 1000 to 2- 104 + 1000, ( Table 


22). 


datal= {{499,162}, {735,129}, {825,121}, {885,114}, {887,113}, 
{892, 108}, {895,105}, {894, 106}, {900, 100}, {903, 97}, 
{897, 103}, {891,99}, {904, 96}, {898, 102}, (908, 92}, 
{902, 98}, {912, 88}, {902, 98}, {908,92}, {902, 98}, 
{906,94}, {902,98}, {912, 88}, {908,92}, {905, 95}, 
{912, 88}, {902, 98}, {908, 92}, {902, 98}, {906, 94}, 
{902, 98}, {912, 88}, {908,92}, {905,95}, {912, 88}, 
{908,92}, {914, 86}, {900, 90}, {908,92}, {912, 88}, 
{925, 75}, {932, 68}, {936, 64}, {930, 70}, {937, 63}, 
{934, 66}, {945, 55}, {925,75}, {940, 60}, {939, 61}, 
{930, 70}, {938, 62}, {946,54}, {949,51} {946,54}, 
{945, 55}, {955, 45}, {944, 56}, {942, 58}, {951,49}, 
{956, 44}, {962, 38}, {953, 47}, {956, 44}, {955, 45}, 
{961, 39} {953,47}, {960, 40}, {955, 45}, {966, 34}, 
{963, 37}, {961,39} {971, 29}, {974, 26}, {963, 37}, 
{965, 35}, {962,38}, {965,35} {970,30} {966, 34}, 
{964, 36}, {970, 30}, {976, 24}, {973,27} {971,29}; 


Table 22: The data.1 


The diagram given below (Fig.13) shows the arrangement of the points at the level (x,y) according to the 
data1[19]. 


Figure 13: The depiction of the intervals False and the primes ones (file:datal) contained in the intervals 
N+06—-—N = 1000. 


Even more macroscopically (Fig.14), the points are shown by the line y = 988.709 — 0.988372 * x, which, 
as they appear, are stacked on its lower right, 

It is obvious that for {x = 999, y = 1.32} and for {y = 1000, « = 0.33}, a value located below the unit and 
means that it can in such a interval, and we are talking about this is for high order integers intervals so 
there is not one prime within the interval 6 = 1000, chosen at random. 


Basic thinking usually in statistics is to calculate at constant shows the number of data we need respec- 
tively, finding the frequencies that appear. This method is a real percentage of realistic Espigue isolates the 
properties showing a group of numbers. 


11.2 The second statistic refers to the count of the primes ones that located in the intervals [N, N + 4] 
with N = 1000 up to 2-104, range 6 = 1000 and the count of the primes within successive ZetaZeros. 
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Figure 14: Graphic depiction of the line y = 988.709 — 0.988372 * « with the archive(datal) points. 


After we found the roots kf of general Equation 6+ Zetazero (kin) -Zetazero (kf) = 0 of given initial 
value kin, usually by the Newton method. The data 19 (Table 23), that we have met with the above method 
are. 


data= {{1000, 162}, {2000, 129}, {3000, 121}, {4000, 114}, {5000, 113}, {6000, 108}, 
{8000, 105}, {9000, 106}, {10000, 106}, {11000,100}, {12000, 97} 
{13000, 103},{14000, 99}, {15000, 96}, {16000, 102}, {17000, 92}, 
{18000, 98}, {19000, 88}, {20000, 98}, {21000, 92},{22000, 98}, 
{23000, 94}, {24000, 98}, {25000, 88}, {26000, 92}, {27000, 95}, 
{27000, 95}, {28000, 88}, {29000, 92}, {30000, 86}, {40000, 90}, 
{50000, 92}, {60000, 88}, {1001000, 75}, {2001000, 68}, {3001000, 64}, 
{4001000, 70}, (5001000, 63}, {6001000, 66}, {7001000, 55}, {8001000, 75}, 
{9001000, 60}, {10001000, 61}, {20001000, 70}, 40001000, 62}, 
{60001000, 54}, {80001000, 51}, {100001000, 54}, {200001000, 55}, 
{600001000, 56},{800001000, 58}, {1000001000, 49}, {10000001000, 44}, 
{20000001000, 38},{30000001000, 49},{10000001000, 44},{20000001000, 38}, 
{30000001000, 47}, {40000001000, 44},{60000001000, 45},{80000001000, 39}, 
{100000001000, 47}, {200000001000, 40},{400000001000, 45},{600000001000, 34}, 
{1000000001000, 37}, {2000000001000, 39},{3000000001000, 29}, 
{4000000001000, 26},{5000000001000, 39}, {4000000001000, 26}, 
{5000000001000, 37}, (6000000001000, 35}, {7000000001000, 38}, 
{8000000001000, 35}, {9000000001000, 30}, {10000000001000, 34}, 
{20000000001000, 36}, {40000000001000, 36}, {70000000001000, 30}, 
{90000000001000, 24}, {110000000001000, 27}, {200000000001000, 29}}; 


Table 23: The data 
Using the NonlinearModelFit[19] method of the function y = (a+ cx*2x)/(log[d+b*a]) and after specifying 
the variables {a, b, c,d} the following result the function will be: 


_ 1108.254246288494 — 1.116325646187134 x 10-V ez 
~ log[—4503.90336177023 + 5.4256351714030812| 


with very good approach and value performance in each pairing.The diagram (Fig.15) of the above equation 
is shown in more detail below with gravity in the latest data. 

y = (a+cxx)/(log[d+b*z]). Thus a test value for r = 2-10 gives us y = 25.6, close to 29 we took with 
the analysis. They follow statistics ANOVA and t-Statistic [19,20] and we get the (Tables 24, 25). 


Analysis of variance (ANOVA) is an analysis tool used in statistics that splits an observed aggregate vari- 
ability found inside a data set into two parts: systematic factors and random factors. The systematic factors 
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Figure 15: The graphic depiction of the function 


have a statistical influence on the given data set, while the random factors do not.An ANOVA test is a way 
to find out if survey or experiment results are significant. In other words, they help you to figure out if you 
need to reject the null hypothesis or accept the alternate hypothesis. 


Analysis-ANOVA 


ss MS 
441724. 110431. 
1329.65 17.9683 
443054. 

71143.8 


Table 24: 


Analysis-T-statisrtic 
Estimate Standard Error t- Statistic P- Value 
1108.25 8.79256 126.044 4.12076 10-** 


-1.11633* 107" 5.63966 107"? -1.97942 0.0514887 
5.42564 0.186181 29.1417 2.68171 x 102 
- 4503.9 0.153555 -29330.8 3.51117 10-2 


Table 25: 


As we can see from the results, we have an important Standard Error only for the variable a. The other 
variables are observed to have a low statistical error and have good compatibility. By adapting the method 
aS we see, we associate two lists of results of the number of primes and long intervals, which although 
disproportionately, work together with impeccable and good contact. 


Ronald Fisher introduced the term variance and proposed its formal analysis in a 1918 article The Corre- 
lation 

Between Relatives on the Supposition of Mendelian Inheritance. His first application of the analysis of vari- 
ance was published in 1921. Analysis of variance became widely known after being included in Fisher’s 1925 
book Statistical Methods for Research Workers. 


One of the attributes of ANOVA that ensured its early popularity was computational elegance. The structure 
of the additive model allows solution for the additive coefficients by simple algebra rather than by matrix 
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calculations. In the era of mechanical calculators this simplicity was critical. The determination of statis- 
tical significance also required access to tables of the F function which were supplied by early statistics texts. 


12. Standard—Error and Confidence—Interval 


The standard error (SE) of a statistic (usually an estimate of a parameter) is the standard deviation of 
its sampling distribution or an estimate of that standard deviation. If the statistic is the sample mean, it 
is called the standard error of the mean (SEM). The sampling distribution of a population mean is 
generated by repeated sampling and recording of the means obtained. This forms a distribution of different 
means, and this distribution has its own mean and variance. Mathematically, the variance of the sampling 
distribution obtained is equal to the variance of the population divided by the sample size. This is because 
as the sample size increases, sample means cluster more closely around the population mean. Therefore, 
the relationship between the standard error of the mean and the standard deviation is such that, for a given 
sample size, the standard error of the mean equals the standard deviation divided by the square root of the 
sample size. In other words, the standard error of the mean is a measure of the dispersion of sample means 
around the population mean. Where we observed, in addition to the initial measurements #11-#12, the 
Standard - Error also the confidence interval stabilized at good and acceptable values with central average 
values around the lateral expected values. The values are given in two consecutive tables that we add below 
and refer to the whole values of range (Table 26) 


served Predicted Standard Error Confidence Interva 


Predicted Standard Error Confidence interva 


62 162.352 4.17901 (154.025, 170.679) 50.5969 0.360493 {49.8786 51-3152} 
129 126.574 0.851191 (124.878, 128.27} 49.941 0.356262 (49.2311. 50.6509} 
22 118.232 0.759511 116.719, 119.745) aDA430 O.253002 (48-7403, 50.1473} 
as 113.637 0.729151 (112.184, 115.09) S72 pagers C419 45, AS AGO 

13 110.529 0.711538 109.112, 111.947) reper eas aos pean 
os 108.21 0.699052 (106.817, 109.602) shat paesseanpeal gai nas asec 
os 104.861 0.68137 (103.503, 106.219} caw cuiec GGacass fa a aes 

06 103.582 0.67462 (102.237, 104.926} aoa ‘gondaacs (40.7576. 41.9523} 
106 102.476 0.668765 (101.143, 103.808) “24> O 207474 (40.4199. 41.6053} 
100 101.504 0.6636 (100.182, 102.826) 39.9228 0.290385 (39.4042, 40.5614} 

100.64 0.6s8982 (99.3265, 101.953} 328.9993 0.223419 (32.4345. 39.564} 
103 99.862 0.65481 98.5573, 101.167) 38.4427 0.279364 {37-226 32.9993} 
99.1565 0.651008 97.8593, 100.454) 37-7577 0.274288 {37-2112, 38.3043} 
98.5114 0.647517 97.2212, 99.8016) 36.8426 0.267924 {36.3147_. 37-3824} 
202 97.9179 0.644293 (96.6341, 99.2017} 36.3208 0.265411 {$5-7919. 36.8496} 
97.3687 0.6413 96.0909, 98.6465) 35.9441 0.26516 {35-4158 36.4725} 
96.858 0.638507 95.5858, 98.1303} amps ae 5 Se ailula aaeaimgeriarnste J 
96.3812 0.635891 (95.1142, 97.6482) soo" SS 5 Spe coder eae 
95.9342 0.633431 94.6721, 97.1964} ae parses eee seo 
95.5138 0.631112 (94.2563, 96.7713) See eas eee eTaS nari 
95.1171 0.628918 93.864, 96.3703) Siegen a (34.1064. 35.2754} 
94.7418 0.626837 93.4928, 95.9908) Saunaa. Riansoos (32.7926. 34.4043} 
94.3858 0.624859 (93.1408, 95.6309) 2227196 0.695498 {30.2338. 33.6054} 

s 94.0474 0.622974 92.8061, 95.2887) 30.6249 1.16905 (28.3555. 33.0143} 
is2 93.7249 0.621176 92.4872, 94.9627} 29.7962 1.42221 {26.8314, 32.7621} 
9s 93.4171 0.619455 92.1829, 94.6514) 28.9642 1.80737 {25.3635_. 32.5661} 

s 93.1228 0.617807 91.8918, 94.3538) 25.s626 3.23502 {19.1167. 32.0025} 
92 92.8408 0.616226 91.613, 94.0687) 

é 92.5703 0.614706 (91.3455, 93.7951) 

90 90.3458 0.602126 89.1461, 91.5456) 
92 88.7017 0.592736 (87.5206, 89.8828) 

8 87.4066 0.585287 86.2403, 88.5728) 

s 71.4688 0.490798 70.4909, 72.4468) 

s 68.4113 0.47218 67.4704, 69.3521) 

4 66.7409 0.461946 65.8205, 67.6614) 

° 65.6044 0.454957 (64.6979, 66.5109) 

3 64.7492 0.449683 (63.8531, 65.6452) 

é 64.0667 0.445467 (63.1791, 64.9544) 

s 63.5009 0.441964 (62.6202, 64.3815) 

s 63.0187 0.438976 62.1441, 63.8934) 

° 62.5995 0.436374 (61.73, 63.469) 

1 62.2291 0.434073 61.3642, 63.0941) 

° S9.8979 0.419533 59.062, 60.7338} 

2 $7.735 0.405955 (S6.9261, 58.5439) 

4 56.5407 0.398419 55.7468, 57.3346) 

1 $S.7229 0.393243 54.9393, 56.5064) 

4 55.1046 0.389322 54.3289, 55.8804) 

s 53.2687 0.377631 (52.5163, 54.0212) 

s 51.5512 0.366632 50.8207, 52.2817) 

Table 26: 
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13. Statistical comparison 


With statistical comparison it appears that as long as we moving to a higher order of integer-size that 
are within a given interval, we have chosen 6 = 1000, the count of the primes diminishes to disappear or to 
there are 1 — 2 primes at high order intervals of more than 107°. This also agrees with of Gram’s law. In 
particular, the problem of distribution of the differences t,41 —t, (that is of difference of ZetaZeros) 
is considered [21]. If we accept Gram’s law then the order of this difference does not exceed the quantity 


thoi — tn & ith) + 0, — oo for much larger integers [20,21] then mean and their mean value is close to, 


’ ee In In(n) 
Fe eS Ta 


> 0,n > oo 


And as the above analysis of-False intervals- as we have shown, it is compatible with this result of 
law Gram’s. 


14. From above analysis they arise 3 big conclusions: 


1%. The number(Count) of primes located within 6 intervals is gradually decreasing with a higher or- 
der size n of 10”. 

2°¢. The distribution of the number(Count) of the primes, within interval 6 follows a Nonlinear correlation 
by the function y = (a+ b- «)/(log(d+c- x) similar to the number (Count) of the primes m(x) = x/log(z), 
which is apply approximating for large numbers. 

3°¢. The measurement(Count) of false intervals and the count of the primes follows a linear correlation and 
it increases if the order of magnitude size of the integer increases. 


Epilogue 


This analysis has proved the Riemann Hypothesis is correct, since the real part of the non-trivial zero- 
equation functions is always constant and equals 1/2. This proven by 3 independent methods. Theorems 
{1.Theorem I, II (pages 1 — 7), 2. Theorem (pages 7 — 8), 3. Theorem (pages 8 — 10)}. This analysis has also 
demonstrated that the imaginary part of the non trivial zero zeta function accepts certain values according 
to the intervals defined by the solution of the 3 based functional equations {Eq.set (I, II), III (page 10)}. 


The Riemann Zeta function as defined by Riemann is ¢(s,0) = ¢(s) and applies Riemann’s functional 
equations. If the equations are solved, the only solution for the real part of s is 1/2. Many other scientists 
refer to other forms such as ¢(s,a) or ¢(Q(x),a) or Dirichlet series and many other cases, which are not 
related to the original definition and therefore pertain to the analysis of other cases. 


<The question is not simply whether Riemann’s Hypothesis is valid or not, but which ¢-functions do we 
accept that the Hypothesis as defined by Riemann includes; should we accept the general view of existing 
¢-functions that the Hypothesis includes or not? So if we take as a general assumption of the Hypothesis all 
parallel ¢-functions then it is not 100% true. Of course it is also true that the Riemann Hypothesis has mul- 
tiple connections with other sciences and gives directions appropriate and grouped according to the function 
that we accept to be valid in each particular application case. If we accept generalized ¢-functions, this is 
taken as the only real answer to the Riemann Hypothesis conjecture that has occupied the mathematical 
community for more than 200 years because it generalizes the hypothesis> 


69 


References 


A Proof of the Riemann Hypothesis, Murad Ahmed Abu Amr (2016, arXiv.org) 
A Geometric Proof of Riemann Hypothesis, Kaida Shi p.20-45,(2003 arXiv.org) 


An Exploration of the Riemann Zeta Function and its Applications p.15-60, 2006. 


Themis Mitsis,(Complex analysis I) p.17-30, 


Mantzakouras Nikos Solve Polynomials and Transcendental equations using Mathematica, 
p.50-80,2012, also in internet a part 


Nikolentzoc Poluqronic, R-H p.9-22, 2008. 


A theory for the zeros of Riemann and other 2014, L-functions Guilherme Frand Andre Le Clairy 
p.38-50 


Titchmarsh E C. The zeros of the Riemann zeta-function. Proc Roy Soc Ser, 1935,151: 234-255. 
Titchmarsh E C. The zeros of the Riemann zeta-function. Proc Roy Soc Ser, 1936,157: 261-263. 


Van de Lune J, te Riele H J J, Winter D T. On the zeros of the Riemann zeta-functionin the critical 
strip. Math. Com, 1986, 46: 667-681. 


Hardy G H, Littlewood J E. The zeros of Riemann’s zeta-function on the critical line. Math Z, 1921, 
10: 283-317. 


H.M. Edwards, Riemann’s Zeta Function, Academic Press, New York, London, 1974. 

D. Zagier, The first 50 million prime numbers”, Mathematical Intelligencer (1977), p.7-19. 
The Mathematical Intelligencer (1977), p.7-19. 

Andrew Odlyzko. The 10?! -st zeros of the Riemann zeta function, 1988. Tables of Zeros. 


Cloizeaux and Mehta, The Widom—Dyson constant for the gap probability in random matrix theory, 
2007, p.26-47. 


Hugh Montgomery, The pair correlation of the zeros function, 1970. 
The “datal & data & NonlinearModelFit Statistic”, have been done with the program mathematica. 


’The 1013 first zeros of the Riemann Zeta function, and zeros computation at very large heigh’, 
Xavier Gourdon, version 2004. 


On the Gram’s Law in the Theory of Riemann Zeta Function. M.A. Korolev, 2010, p.4-6. 


Euler Gamma function page 123, http://www.math.leidenuniv.nl/~evertse/ant13-8. pdf 
Video of Enrico Bombieri, https: //www.ias.edu/video/analysis/davenport-hellbron 


Solve Polynomial and transcendental Equations with use Generalized Theorem (Method Lagrange), 


arxiv.org/search/?query=mantzakourast+&searchtype=all&source=header), 
Mantzakouras Nikos 2021 


70 


25 


26 
27 
28 


29 


30 


31 


32 


33 


Zeros of Lattice Sums: 1. Zeros of the Critical Line R.C. McPhedran, 
https://arxiv.org/pdf/1601.01724. pdf 


The inverse Riemann zeta function, A Kawalec (p. 61-63).18 June 2021. 
E. Saias, A. Weingartner, Zeros of Dirichlet series with periodic coefficients, Acta Arith. 140 (2009) 


On the concavity of Dirichlet’s eta function and related functional inequalities.Horst Alzer,Man 
KamKwong,2015 


Hafner, J.: Zeros on the critical line for Dirichlet series attached to certain cusp forms. Math. Ann. 
264, 2t-37 (1983) 


Heath-Brown, D.R.: Simple zeros of the Riemann zeta-function on the critical line. Bull. London 
Math. Soc. 11, 17-18 (1979) 


Levinson, N.: More than one-third of zeros of Riemann’s zeta-function are on t = 1/2. Adv. Math. 
13, 383-436 (1974) 


Davenport, H.; Heilbronn, H., On the zeros of certain Dirichlet series I, I. J. London Math. Soc. 11 
(1936), 181-185, 307-312. 


Spira, Robert., Some zeros of the Titchmarsh counterexample. Math. Comp. 63 (1994), no.208, 
747-748. MR1254148 (95a:11081) 


71 


